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ABSTRACT 

Recent work suggests that the mass-loss rate of the primary star t] A in the massive colliding 
wind binary rj Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results 
from large- (±1545 au) and small- (±155 au) domain, 3D smoothed particle hydrodynamics 
(SPH) simulations of i] Car’s colliding winds for three /; A mass-loss rates ( /V/, (A = 2.4, 4.8 
and 8.5 x 1 0 4 Mq yr 1 ), investigating the effects on the dynamics of the binary wind-wind 
collision (WWC). These simulations include orbital motion, optically thin radiative cooling 
and radiative forces. We find that M, u greatly affects the time-dependent hydrodynamics 
at all spatial scales investigated. The simulations also show that the post-shock wind of the 
companion star tjb switches from the adiabatic to the radiative-cooling regime during periastron 
passage (tp ~ 0.985-1.02). This switchover starts later and ends earlier the lower the value of 
M, u and is caused by the encroachment of the wind of rj A into the acceleration zone of jjb’s 
wind, plus radiative inhibition of ?) B ’s wind by ri A . The SPH simulations together with ID 
radiative transfer models of rj A s spectra reveal that a factor of 2 or more drop in M, lA should 
lead to substantial changes in numerous multiwavelength observables. Recent observations are 
not fully consistent with the model predictions, indicating that any drop in M , ]A was likely by 
a factor of <2 and occurred after 2004. We speculate that most of the recent observed changes 
in tj Car are due to a small increase in the WWC opening angle that produces significant 
effects because our line of sight to the system lies close to the dense walls of the WWC zone. 
A modest decrease in M nA may be responsible, but changes in the wind/stellar parameters of 
0b, while less likely, cannot yet be fully ruled out. We suggest observations during i] Car’s 
next periastron in 2014 to further test for decreases in M, u . If M, lA is declining and continues 
to do so, the 2014 X-ray minimum should be even shorter than that of 2009. 
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ing wind binary with a highly eccentric (e ~ 0.9), 5.54 yr orbit 
(Damineli, Conti & Lopes 1997; Hillier et al. 2001: Damineli et al. 
2008a, b; Corcoran et al. 2010). The primary of the system, p A , is 
our closest (2.3 ±0.1 kpc, Smith 2006) example of a very mas- 
sive star. A luminous blue variable (LBV) with a current mass 
of ~T00Mq and a powerful stellar wind with a mass-loss rate 
M ^ 8.5 x 10~ 4 MQyr~' and terminal speed Uoo 420kms~' 
(Hillier et al. 2001, 2006; Groh et al. 2012a, hereafter H01, H06 
and G12a, respectively), ij A is in a relatively short and poorly 
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nae- stars: winds, outflows. 

1 INTRODUCTION 

Deep within Eta Carinae’s spectacular bipolar ‘Homunculus’ neb- 
ula lies an extremely luminous (Lxotai > 5 x 10 6 Lq) collid- 
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understood stage of stellar evolution. Such very massive stars can 
lose substantial amounts of mass in giant, violent outbursts as they 
rapidly transition from hydrogen burning to helium core burning 
(Conti 1984; Humphreys & Davidson 1994). The most prominent 
illustration of this is the mid- 1840s ‘Great Eruption’ of Eta Carinae 
(Xi Car), wherein ~10-20 Mq was ejected over a period of ~5— 
15 yr, forming the Homunculus nebula (Davidson & Humphreys 
1997; Smith 2006). Evidence indicates that such strong mass-loss 
dominates the evolution of the most massive (>40 Mq) stars (see 
e.g. Smith & Owocki 2006). The physical mechanisms leading to 
these eruptions are not well understood, in either single stars (Smith 
& Owocki 2006; Smith & Townsend 2007; Quataert & Shiode 2012) 
or binaries (Ofek et al. 2007; Smith 201 1). Determining their cause 
is of fundamental importance for understanding how very massive 
stars evolve into supernovae (SNe), gamma-ray bursts and black 
holes. 

Ground- and space-based multiwavelength observations over the 
last two decades show that rj A has a less-luminous binary compan- 
ion, ? 7 b (Damineli et al. 1997; Ishibashi et al. 1999; Smith et al. 
2004b; Whitelock et al. 2004; Damineli et al. 2008a, b; Gull et al. 
2009; Fernandez-Lajus et al. 2010; Corcoran et al. 2010, here- 
after CIO). Some of the strongest evidence for binarity comes from 
extended X-ray monitoring of the system (Ishibashi et al. 1999; 
Corcoran 2005; Hamaguchi et al. 2007; CIO), r) Car’s Rossi X-ray 
Timing Explorer ( RXTE ) light curve is characteristic of a highly 
eccentric colliding wind binary, the variable X-ray emission arising 
in a wind-wind collision (WWC) shock zone formed between the 
stars (Pittard & Corcoran 2002; Parkin et al. 2009, 2011, hereafter 
P09 and PI 1; CIO). The hard X-ray spectrum ( kT ^4— 5keV) 
requires that rj B have an unusually high wind terminal speed of 
~3000kms -1 , while modelling of the X-ray spectrum suggests an 
enormous mass-loss rate of ~10 -5 Mgyr -1 (Pittard & Corcoran 
2002; P09). Due to p A ’s immense brightness, has never been 
directly observed. Nonetheless, studies of the photoionization ef- 
fects of r ) b on the circumstellar ejecta known as the ‘Weigelt blobs’ 
(Weigelt & Ebersberger 1986) have helped constrain ijg’s temper- 
ature (T e( f ~ 37 000-43 000 K) and luminosity (logi/Lg ~ 5-6, 
Verner, Bruhweiler & Gull 2005; Mehner et al. 2010). 

Monitoring of r) Car with RXTE and Swift shows that, broadly 
speaking, the X-ray flux repeats itself from orbit to orbit with more 
consistency than is shown in any other waveband (CIO; Corco- 
ran 2011). Surprisingly, the observed X-ray minimum in 2009 was 
~50 percent shorter than the minima of 1998 and 2003.5 (fig. 3 of 
CIO). The sudden recovery of the X-ray flux in 2009 probably rep- 
resents a change in the cooling of the WWC shock due to currently 
unknown reasons (P09; Pll). Simple calculations suggest a large 
reduction in p A ’s M, of a factor of ~2-4, is needed to match the 
reduced 2009 minimum (CIO). 

Spectroscopic changes observed in stellar- wind emission features 
such as Ha are also interpreted as being due to a secular decrease 
of t) a ’s mass-loss rate (by a factor of 2-3 between 1999 and 2010, 
Mehner et al. 2010, 2011, 2012). However, monitoring of Balmer 
H<5 wind lines shows no significant changes in the inner wind of rj A 
(Teodoro et al. 2012). Moreover, while a secular change is visible 
in various wind profiles for a direct view of the central source, little 
to no change is seen in wind profiles at high stellar latitudes or 
reflected off of the Weigelt blobs (Gull et al. 2009; Mehner et al. 
2012). Model stellar spectra also predict that a factor of 2 or more 
change in the M of ;? A should lead to very significant changes in 
the observed spectrum (H06). Hei wind lines should become very 
bright, while lines of Fe ii should nearly disappear (fig. 1 1 of H06). 


Such changes have not been observed, and so the puzzle of what 
recently occurred in Car remains. 

Here, we present results from large- (±1545au) and small- 
(±155 au) domain, 3D smoothed particle hydrodynamics (SPH) 
simulations of i] Car’s massive colliding winds for three different 
mass-loss rates of ;y A (2.4, 4.8 and 8.5 x 10~ 4 Mg yr -1 ). The goal 
is to investigate how decreases in ?y A ’s M affect the 3D geometry 
and dynamics of ?; A ’s optically thick wind and spatially extended 
WWC regions, both of which are known sources of observed X-ray, 
optical, UV and near-infrared (IR) emission and absorption (Pit- 
tard & Corcoran 2002; van Boekel et al. 2003; Martin et al. 2006b; 
Nielsen et al. 2007; Weigelt et al. 2007; Gull et al. 2009, 2011; 
P09, CIO; Groh et al. 2010a, b; Pll; G12a; Madura & Groh 2012; 
Madura et al. 2012; Teodoro et al. 2012). We use several compu- 
tational domain sizes in order to better understand how the various 
observables that form at different length-scales may be influenced 
by declines in ;; A ’s M. This is the first such parameter study of 
r i Car in 3D, which is essential since orbital motion greatly affects 
the shape and dynamics of the WWC zones during periastron pas- 
sage (Okazaki et al. 2008, P09; Pll; Madura & Groh 2012; Madura 
et al. 2012). The simulations in this paper can help constrain r/ Car’s 
recent mass-loss history and possible future state. 

The following section describes the 3D SPH simulations, while 
Section 3 presents the results. A qualitative discussion of the ob- 
servational implications is in Section 4. Section 5 summarizes our 
conclusions and outlines the direction of future work. 

2 THE 3D SPH SIMULATIONS 

The hydrodynamical simulations in this paper were performed with 
an improved version of the 3D SPH code used in Okazaki et al. 
(2008) and Madura et al. (2012, hereafter M12), which is based 
on a version originally developed by Benz et al. (1990) and Bate, 
Bonnell & Price (1995), to which we refer the reader for further 
details. The linear and non-linear SPH artificial viscosity parame- 
ters are «sph = 1 and jSsph = 2, respectively (Okazaki et al. 2008). 
Optically thin radiative cooling is now implemented using the Ex- 
act Integration scheme of Townsend (2009). The radiative cooling 
function A(7) is calculated with cloudy 90.01 (Ferland et al. 1998) 
for an optically thin plasma with solar abundances. The pre-shock 
stellar winds and rapidly cooling dense gas in the WWC region are 
assumed to be maintained at a floor temperature = 10 4 K due to 
photoionization heating by the stars (Pll). The same initial tem- 
perature is assumed for both winds for simplicity. The effect of the 
initial wind temperature (r w i nt i) on the flow dynamics is negligible 
(Okazaki et al. 2008). 

Radiative forces are incorporated in the SPH code via an ‘anti- 
gravity’ formalism developed by two of us (ATO and SPO), the 
details of which can be found in Appendix A and Russell (2013). 
We parametrize the individual stellar winds using the standard ‘beta- 
velocity law" v(r) — 1 ) 00(1 — R*/r)^, where Uoo is the wind terminal 
velocity, R „ the stellar radius and f) a free parameter describing the 
steepness of the velocity law. Effects due to ‘radiative braking,’ in 
which one of the stellar winds experiences a sudden deceleration 
before reaching the WWC zone (Gayley, Owocki & Cranmer 1997, 
Pll) are not included. Such braking effects are not expected to play 
a prominent role in r] Car (P09; Pll; Russell 2013). We include 
the more important velocity-altering effects of ‘radiative inhibition’ 
(RI), in which one star’s radiation field reduces the net rate of ac- 
celeration of the opposing star’s wind (Stevens & Pollock 1994, 
P09; Pll). However, because we fix the stellar mass-loss rates in 
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our antigravity approach, possible changes to the mass-loss due to 
RI are not included. These mass-loss changes are not expected to 
be significant in q Car and should not greatly affect our results or 
conclusions (see Appendix A3 for details). 

Using a standard xyz Cartesian coordinate system, the binary orbit 
is set in the xy plane, with the origin at the system centre-of-mass 
(COM) and the major axis along the x-axis. The two stars orbit 
counter-clockwise when viewed from the ±j;-axis. By convention, 
1 = 0 (4> = t/ 2024 = 0) is defined as periastron. Simulations are 
started with the stars at apastron and run for multiple consecutive 
orbits. Orbits are numbered such that (j> = 1.5, 2.5 and 3.5 corre- 
spond to apastron at the end of the second, third and fourth full 
orbits, respectively. 

The outer spherical simulation boundary is set at r = 10a and 
100a from the system COM for the small- and large-domain simu- 
lations, respectively, where a = 15.45 au is the length of the orbital 
semimajor axis. Particles crossing this boundary are removed from 
the simulations. The total number of SPH particles used in the 
small- and large-domain simulations is roughly the same (between 
~5 x 10 5 and 9 x 10 5 , depending on the value of the primary 
M). The small-domain simulations are better suited for studying 
the complex physics and geometry of the WWC zone very close 
to the central stars where most of the instantaneous changes ob- 
served in numerous spectral lines, the optical/near-IR/radio contin- 
uum fluxes and X-rays across q Car’s 5.54 yr orbital cycle occur 
(see e.g. Damineli et al. 2008b). 

The large-domain simulations are comparable in size to past and 
planned Hubble Space Telescope/Space Telescope Imaging Spec- 
trograph ( HST/STIS ) mapping observations of the interacting stellar 
winds in q Car’s central core (~ ± 0.67 arcsec Rt ±1540 au, Gull 
et al. 2011; M12; Teodoro et al. 2013). As demonstrated by Gull 
et al. (2011) and M12, 3D simulations at this scale are necessary for 
understanding and modelling the extended, time-variable forbidden 
line emission structures that are spatially and spectrally resolved by 
HST/ STIS. We note that the SPH formalism is ideally suited for such 
large-scale 3D simulations due to the significantly decreased com- 
putational cost, as compared to grid-based hydrodynamics codes 
(Price 2004; Monaghan 2005). 

The adopted simulation parameters (Table 1 ) are consistent with 
those derived from the available observations, although there has 
been some debate on the exact value of q A ’s mass-loss rate (here- 
after M va ). For the current study, we investigate three values of 
M, a , which we refer to as Cases A— C. The largest adopted value, 
8.5 x 10~ 4 Mq yr~' (Case A), is based on the latest cmfgen radia- 
tive transfer model fits to HSTI STIS spatially resolved spectroscopic 
observations of q Car (G12a). H01, H06, Davidson et al. (1995) and 
Cox et al. (1995) derived a mass-loss rate of 10~ 3 Mq yr~' and a 
wind terminal speed of 500 km s” 1 for q A . The good match be- 
tween the synthetic [Fern] spectroimages of M12 and HST/ST1S 
observations taken between 1998 and 2004 (Gull et al. 2009) fur- 
ther suggests M 1]A was ~10~ 3 MQyr _1 during that time. G12a 
improved upon the original model fits of H01, H06 by using a 
slightly lower M, 1A (8.5 x 10~ 4 MQyr~') and Uoo (420km s _1 ). 
The optical and UV spectra modelled by G12a were taken by HST 
in 2000 March and 2001 April, corresponding to a phase near apas- 
tron. We adopt the G12a value and assume it as the initial M, 1A 
before any possible drop after 2001. We also use the same for 
Pa (Voo,r, A = 420kms _1 ) in each simulation. 

Our second adopted M, 1A (4.8 x 10~ 4 Mq yr -1 , Case B) is con- 
sistent with that derived by P09 and used in the 3D adaptive mesh 
refinement (AMR) hydrodynamical simulations of PI 1 . The lowest 
value of M, 1a used in our simulations (2.4 x 10~ 4 Mq yr -1 , Case C) 


Table 1 . Stellar, wind and orbital parameters of the 3D SPH 
simulations. 


Parameter 

PA 

VB 

Reference 

M„ (Mq) 

90 

30 

H01, 008 

/?. <Rq) 

60 

30 

H01, H06 

Twind (10 4 K) 

3.5 

3.5 

008, this work 

M ( 10 — 4 Mq yr~ 1 ) 

8.5, 4.8, 2.4 

0.14 

G12a, P09 

Uoo (km s -1 ) 

420 

3000 

G12a, PC02 

P 

1 

1 

H01, G12a 

n 

0.12, 0.21,0.42 

this work 

q 

22.28, 15.22, 9.59 

H01, P09, this work 

Porb (d) 

2024 


D08a 

e 

0.9 


C01.P09 

a (au) 

15.45 


CO 1,008 


Notes. M t and R t are the stellar mass and radius. T wln( i is the ini- 
tial wind temperature. M, Uqo and fS are the stellar-wind mass- 
loss rate, terminal speed and velocity-law index, respectively. 
r ; = (Muoo)ij B /(3Tu 00 ), A is the secondary/primary wind momentum 
ratio, q is the primary/secondary stellar luminosity ratio assuming 
L„ jf A = 5 x 10 6 Lq, P or b is the period, e is the eccentricity and a is 
the length of the orbital semimajor axis. 

References: C01: Corcoran et al. (2001); H01: Hillier et al. (2001); 
PC02: Pittard & Corcoran (2002); H06: Hillier et al. (2006); D08a: 
Damineli et al. (2008a); 008: Okazaki et al. (2008 ); P09: Parkin et al. 
(2009); G12a: Groh et al. (2012a). 

is nearly identical to the M, 1A obtained by Pittard & Corcoran (2002) 
in their attempts, using 2D hydrodynamical simulations of the col- 
liding winds, to fit a Chandra grating spectrum of q Car collected 
near apastron. This same M, A was later used by Okazaki et al. 
(2008) in a 3D SPH simulation to model q Car’s RXTE X-ray light 
curve. We note that all of the X-ray modelling studies mentioned 
assume wind terminal speeds of 500 and 3000 km s _1 for q A and 
q B , respectively, and that we use the same M and as PI 1 for q B 
(Table 1). 

3 RESULTS 

3.1 Small-domain simulations 

We begin with the results of our r = 10a simulations, focusing on 
slices in the orbital plane. Results for slices in the xz and yz planes 
are in Appendix B 1 . Movies of the simulations are also available in 
the online version of this paper (see Supporting Information). 

Figs 1 and 2 display the density and temperature at five select 
orbital phases, demonstrating the key features of the WWC and the 
effects of orbital motion. As in previous 3D simulations (Okazaki 
et al. 2008; PI 1; Madura & Groh 2012; M12), the fast, low-density 
wind of q B carves a large cavity out of the slower, denser wind of q A 
for most of the 5.54 yr orbit. Around apastron (row a of both figures) 
when orbital speeds are their lowest (~ 14 km s~ 1 for q B with respect 
to the system COM), this cavity and the WWC zone maintain the 
expected axisymmetric conical shape. The measured half-opening 
angle of each shock cone (55°, 65° and 75° for q fb 0.12, 0.21 and 
0.42, respectively) is consistent with that predicted by the analytic 
asymptotic formula (28) of Canto, Raga & Wilkin (1996) (53?5, 
62? 4 and 74? 2) and increases as the value of M, 1A decreases. This 
increase is due to the changing wind momentum balance, which 
also moves the apex of the WWC zone closer to q A . At apastron, 
the WWC apex is ~22 au (18 au) from q A when q ~ 0.12 (0.42). 

The WWC zone in each simulation consists of a distended shock 
containing hot (T > 10 6 K), low-density q B wind material and a 
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Case B 

M = 4.8 x 10 4 M 0 /yr 


a 

0 = 0.50 


b 

0 = 0.90 


c 

0 = 1.00 


d 

0 = 1.03 


e 

0 = 1.10 


j'. 





-20 18 -16 -14 

X log density [g/cm 3 ] 


Case C 

M = 2.4 X 10 4 M q /yr 


-10 -5 0 5 10 



Figure 1. Slices in the orbital plane from the small-domain 3D SPH simulations of r] Car (Table 1) at orbital phases 0 = 0.5 (apastron), 0.9, 1.0 (periastron), 
1.03 and 1.10 (rows, top to bottom). Columns correspond to the three assumed ri a mass-loss rates. Colour shows log density in cgs units. The spherical 
computational domain size is r = 10a ^ 155 au ~0.067 arcsec (D = 2.3 kpc). Axis tick marks correspond to an increment of la ^ 15.45 au. The orbital motion 
of the stars is counter-clockwise. is to the left and ?i b is to the right at apastron. 
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Case A 

M = 8.5 x 10 4 M 0 /yr 


Case B 

M = 4.8 x 10 4 M 0 /yr 


Case C 

M = 2.4 x 10 4 M g /yr 


a 

0 = 0.50 


b 

0 = 0.90 


t 




A 


0 = 1.00 


d 

0 = 1.03 


e 

0 = 1.10 




Figure 2. Same as Fig. 1, but with colour showing log temperature. 


much thinner, colder ( T 10 4 K) region of dense post-shock tj a 

wind material, separated by a contact discontinuity (CD). The post- 
shock density of the wind is roughly an order of magnitude 
larger than the pre-shock density. The post-shock density of the ijb 


wind is approximately a factor of 4 higher than the pre-shock den- 
sity, as expected for adiabatic shocks. The slower, denser r) A wind 
radiates more efficiently and has a characteristic cooling parame- 
ter / = tcooi/lescape = Vgdn/M-i 1 (where u 8 is the pre-shock 
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wind velocity in units of 1000 km s~', dn is the distance to the 
CD in units of 10 l: cm and M_ 7 is the mass-loss rate in units of 
10~ 7 Mq yr -1 , see Stevens, Blondin & Pollock 1992). The lower 
density, high-velocity r ) B wind has / > 1 and is adiabatic for most 
of the orbit. The maximum expected post-shock temperature in the 
r/ B wind (T = 3mu 2 /16k, where m = 10~ 24 g is the average mass 
per particle assuming solar abundances, v w is the pre-shock wind 
speed and k is Boltzmann’s constant, Stevens et al. 1992) at apas- 
tron is ~ 1.2 x 10 s K, which is very near what is observed in the 
simulations at the WWC apex (~ 10 s K, Fig. 2). 

The overall fraction of jjb’s wind shocked to high temperatures 
increases with M, IA . The areas where the post-shock r) B wind is 
hottest (in bright yellow in rows a and b of Fig. 2) increase in size 
as M, a is decreased. This is because for a given pre-shock wind 
speed, oblique shocks (like those in Case A) are less efficient at 
thermalizing the flow and produce lower post-shock temperatures 
(Pittard 2009). In contrast, Cases B and C have WWC zones that are 
more normal to the pre-shock flow, resulting in higher post-shock 
temperatures in the outer wings of the WWC region. 

The post-shock ij A wind region appears to become somewhat 
thinner and less dense the lower the value of M 1]A . Since the col- 
liding winds have very different speeds there is a velocity shear 
across the CD that can excite Kelvin-Flelmholtz (KH) instabilities 
(Stevens et al. 1992; PI 1). The thin, dense and rapidly cooling post- 
shock ij a wind is also subject to non-linear thin-shell instabilities 
(NTSIs, Vishniac 1994; PI 1). As M r]A is decreased, the WWC zone 
seems to become more unstable. One possible reason for this is be- 
cause the time-scale for exponential growth of the KH instability is 
oc ^/PiPi/ipi + P2)' where p\ and p 2 are the densities on the two 
sides of the CD (Stevens et al. 1992). Therefore, as M nA is lowered, 
Pi decreases and the growth time-scale for the KH instability short- 
ens. We note, however, that standard SPH schemes are notorious 
for under-resolving the KH instability (Agertz et al. 2007; Price 
2008), so these results should be interpreted with caution. In the 
case of the NTSI, the overall stability of a dense shell depends on 
the shell's thickness (Vishniac 1994; Blondin & Koerwer 1998). 
The thicker, denser compressed ‘shell’ of primary wind in Case A 
is thus probably less prone to the NTSI. 

The aberration angle and degree of downstream curvature of the 
WWC zone are determined by the ratio of the orbital speed to the 
pre-shock wind speed (Parkin & Pittard 2008; Pittard 2009). As 
the stars approach periastron, the orbital speed of r] A relative to )/ B 
(^360 km s~ 1 ) increases to a value near its wind speed (420 km s _1 ), 
thus highly distorting the WWC zone (Fig. 1, rows b and c). The 
increasing orbital speeds furthermore cause the post-shock rj B gas 
in the leading arm of the WWC zone to be heated to higher tem- 
peratures than the gas in the trailing arm (rows b and c of Fig. 2), a 
result also found by PI 1. 

Because of the high orbital eccentricity, at tf) ~ 0.97—1.0, ;y B 
and the WWC apex become deeply embedded in ij A ’s dense wind 
(row c of Figs 1 and 2). The rapid decrease in stellar separation at 
these phases and the adopted f = 1 velocity laws cause the two 
winds to collide before r] B s wind can reach its terminal speed. 
This, combined with RI of r/ B ’ s wind by r] A , leads to a significant 
decrease in the pre-shock velocity of ?? B ’s wind (compared to phases 
around apastron, see Appendix A). The associated drop in r) B ’s wind 
momentum moves the WWC apex closer to r) B and decreases the 
WWC-zone’s opening angle. Because of these effects, the trailing 
wind of i) B is unable to collide with ?? A ’s downstream wind. As a 
result, the post-shock r] B wind in the trailing arm of the WWC zone 
cools below 10 6 K, making the hot shock noticeably asymmetric 
just before and at periastron (Fig. 2, row c). The higher the value of 


M JtA , the earlier this process starts, with the hot gas in the trailing 
arm vanishing at <p^ 0.985 (0.995) inCase A (Case C). There is also 
a ‘collapse’ of the WWC zone around periastron (see Section 3.1.1 
for details). 

Rapid orbital motion during periastron passage distorts the WWC 
region, giving it a distinct spiral shape that is clearly visible at 
</> = 1.03 (Figs 1 and 2, row d). The wind of rj B also carves a 
narrow, low-density spiral cavity in the back side of ?j A ’ s wind. 
This cavity flows in the — x and — y directions following periastron 
(rows d and e of Figs 1 and 2). The width of the cavity increases 
with decreasing M rlA due to the change in the wind momentum 
balance. The high-velocity tj B wind in the cavity also continuously 
ploughs into the denser, slower ;; A wind that flows in the — x and — y 
directions. This makes the compressed wall or ‘shell’ of ?j A wind 
material grow in density and accelerate to speeds a few hundred 
km s _1 greater than v 00 ^ A . The lower the M Va , the greater the 
outflow speed of the shell due to the more similar wind momenta. 
As the shell sweeps up more mass, it gradually slows to velocities 
comparable to Uoo,ij a . The thickness and density of the spiral shell 
of r/A wind also decrease with decreasing M nA . 

At (p ^ 1.03, the dense layer of post-shock r) A wind located in 
the trailing arm of the WWC region is photo-ablated by the intense 
stellar radiation fields. PI 1 originally identified this photo-ablation 
in their 3D AMR simulations. Interestingly, only in Case A does the 
photo-ablation cut off or ‘snip’ the spiral tail of the r) B wind cavity 
(row d of Figs 1 and 2). This snipping off of the cavity prevents the 
r/ B wind from colliding with and heating the downstream r] A wind, 
causing a noticeable asymmetry in the temperature structure of the 
post-shock r] B gas in the spiral cavity. In Case A, only )/ B wind 
that lies in the cavity in the forward direction remains hot, while 
material in the trailing direction cools to T~ 10 4 K. In Cases B and 
C, the entire ;y B wind cavity is hot, with T > 10 6 K. The temperature 
asymmetry seen in Case A remains as the gas flows outward and is 
also visible at <j> = 1.1 (Fig. 2. row e). 

Recovery of the WWC zone following periastron takes longer the 
higher the value of M, A . In Case A, the hot j/ B -wind shock does not 
fully return until <p ss 1.025, whereas in Case C, the hot shock returns 
by 4> ss 1.015 (see Section 3.1.1 for details). After periastron, ij b 
moves in front of t) A , its wind colliding with and heating the denser 
wind of r/A that flows unobstructed in the +x, —y direction (row 
e of Figs 1 and 2). At <f> ^ 1.1, the arms of the WWC zone are 
so distorted by orbital motion that the leading arm collides with 
the old trailing arm formed before periastron passage. The lower 
the value of M rlA , the sooner this collision takes place. The leading 
arm of the WWC zone, including the portion that collides with the 
old trailing arm, helps form another compressed shell of ;? A wind 
that flows in the +x and — y directions after periastron (row e). The 
higher the ii/, A , the greater the mass in this shell and the slower it 
travels outward. Eventually, orbital speeds decrease and the WWC 
zone and ?; B wind cavity regain their axisymmetric conical shape as 
the system moves slowly back towards apastron. 


3.1.1 'Collapse’ of the WWC zone during periastron passage 

P09 and Pll showed that RI likely reduces the pre-shock velocity 
of )j B ’s wind during periastron passage, resulting in lower post- 
shock temperatures and rapid cooling of the post-shock r/ B gas. 
Such rapid cooling can catastrophically disrupt the WWC zone via 
NTSIs, leading to a ‘collapse’, with dense fragments of post-shock 
y]a wind being driven into the acceleration zone of r] B s wind and 
sometimes colliding with tj B (PI 1). It has been postulated that such 
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a WWC collapse could explain 77 Car’s observed extended X-ray 
minima in 1998 and 2003, with the shorter 2009 minimum caused 
by the lack of a collapse for some unknown reason (P09; C 1 0; P 1 1 ; 
Russell 2013). 

3D AMR simulations of 77 Car by PI 1 surprisingly did not show 
the expected collapse of the WWC zone around periastron. Rapid 
orbital motion, which increases the pre-shock wind velocity and 
post-shock pressure of 77 b’s wind, is one key factor that helped 
prevent a collapse (P 11). A second likely factor is the relatively 
weak coupling between )y A ’s radiation field and 773 ’s wind adopted 
for the simulations. PI 1 suggest that different couplings could result 
in stronger inhibition that leads to a disruption and collapse of the 
WWC zone on to tj B . Due to the high computational cost of 3D 
AMR simulations, Pll did not investigate other radiation-wind 
couplings. They were thus unable to demonstrate a WWC collapse 
around periastron in a 3D simulation that includes orbital motion. 

Since orbital motion, radiative cooling and the velocity-altering 
effects of RI are included in our simulations, they provide a means 
to investigate further the possibility of a WWC collapse around 
periastron. Here, we focus on the inner ±1 a ~ ±15.45 au region 
in the orbital plane during periastron passage {ij> = 0.97—1.03). In 
order to ensure that we are adequately resolving the WWC shocks 
and any possible WWC collapse, we performed a series of very 
small domain 3D SPH simulations with a computational domain 
size of r = 1.5a and eight times the number of particles used in the 
r= 10 a simulations, leading to a roughly factor of 2 improvement in 
the overall spatial resolution compared to the r = 10 a simulations. 
We begin with the M^ K = 8.5 x 10 ~ 4 Mg yr~' simulation results 
(Fig. 3). 

At <f> = 0.98 (Fig. 3, row a), we find that the pre-shock speed 
of 773 ’s wind along the line of centres is significantly reduced to 
'~1950kms~ 1 . The stellar separation at cj> = 0.98 (r 5 au) is small 
enough for the two winds to collide while ? 7 B ’s wind is still ac- 
celerating. The expected pre-shock speed of 773 ’ s wind along the 
line of centres at </> = 0.98 when assuming f) — 1 wind-velocity 
laws and taking into account orbital motion is ~2770kms _1 . RI 
of 773 ’s wind by 77 A explains the difference. Using the analysis in 
Appendix A2, the computed expected total 773 wind speed along 
the line of centres (vi, tot) at the location of ram pressure balance 
r b R* 1.03 au (equation A19) is Rl953kms _1 , in excellent agree- 
ment with the 1950 km s~' found in the SPH simulation. Even at 
apastron, we find that 773 ’s pre-shock wind speed in the simulation 
is reduced to ~2850 km s _1 , compared to the expected 2944 km s _1 
when using simple 7 S = 1 velocity laws (Fig. A3). 

Due to the lower wind speed and shorter distance to the CD at 
ip = 0.98, the colliding ij B wind heats to lower temperatures. The 
shorter distance to the CD also increases the pre- and post-shock 
tj B wind densities. The cooling parameter x in the post-shock ij B 
wind near the WWC apex drops to ~1.6 (Fig. 4). The post-shock 
temperature in the 77 B wind at the WWC apex is reduced, which 
decreases the post-shock thermal pressure that helps support the 
WWC zone (Pll). The lowered r] B wind speed also alters the wind 
momentum ratio and decreases the WWC opening angle. 

Between 0 = 0.98 and 0.99, the stellar separation drops to ~3 au 
and the winds collide at progressively deeper locations within the 
acceleration zone of ? 7 B ’s wind. RI of 773 ’s wind by 77 A also becomes 
stronger (Fig. A3). During this time, the pre-shock 773 wind speed 
along the line of centres in the simulation drops to ~1500kms _1 
(Fig. 3). The WWC opening angle decreases further, the WWC 
apex moves closer to 773 , / drops to well below unity (Fig. 4) and 
radiative cooling takes over. The post-shock 773 wind near the WWC 
apex now cools quickly to 7^ 10 4 K. At this point, we observe 


a ‘collapse’ (P09; Pll) or ‘discombobulation’ (Davidson 2002; 
Martin et al. 2006b) of the inner WWC region. Since a WWC zone 
is still present (although highly unstable due to the two radiative 
shocks), we believe it is more appropriate to say that the post-shock 
ri B wind undergoes a ‘cooling-transition" phase. The entire apex 
and old trailing arm of the hot tj B -shock vanish (Fig. 3, rows b 
and c) as the post-shock rj B gas switches from the adiabatic to the 
radiative-cooling regime. With the thermal pressure from the hot 
ri B wind gone, the cold, dense 77 A wind can penetrate deep into 773 ’s 
wind and come very close to the star, ~3.5 R* , JB (~0.5au) from 
77 b ’s surface at <j) = 0.99. 

At periastron, there is almost a complete collapse of the apex of 
the WWC region in the simulation (row d of Fig. 3). The dense 
wind of 77 A directly between the stars comes within ~1.5 of 
ri B ’s surface. 17 B can only effectively drive a wind in directions away 
from 77 a- Thus, only regions downstream where 773 ’s wind collides 
with the highly distorted spiral leading arm of the dense WWC 
zone are shock-heated to high T (yellow areas in second column of 
Fig. 3, row d). However, these regions have very low densities and 
volumes and are enshrouded by dense primary wind. 

The equations in Appendix A2 predict that r) B should not be 
able to drive a wind along the line of centres at periastron due to the 
strong RI by r] A . In our SPH simulations, however, we launch the 773 
wind at a value of twice the local sound speed (see Appendix Al). 
This, combined with the large transverse velocity of the stars at 
periastron, produces enough of an 773 wind to create a ram pressure 
balance at periastron in our simulations. In the absence of other 
forces (e.g. radiative braking), a full WWC collapse should occur 
along the line of centres at periastron. 

The disappearance of the hot, inner WWC zone (i.e. the total 
cooling-transition phase) lasts ~2.5 months (4> R* 0.983 — 1.022) in 
the simulation of Fig. 3. While the 773 wind speed at (f) = 1.01 is 
enough to restore a WWC at a few R t above 773 ’s surface (row e of 
Fig. 3), it is not until <t> R 1 .02 that 773 moves far enough from r/ A that 
the hot, inner conical shock can start to be fully restored (Fig. 3, row 
f). Therefore, the hottest plasma near the WWC apex responsible for 
the observed X-ray emission is absent from 0 R* 0.983— 1.022. The 
pre-shock speed of i] A ’s wind along the line of centres is also slightly 
reduced by ~30— 40kms~' just after periastron (<fi ~ 1.01 — 1.03) 
due to the rapid orbital motion. 

Figs 4 and 5 show that even with lowered r] A mass-loss rates, 
a cooling-transition phase occurs during periastron passage in the 
post-shock 77 b wind located near the WWC apex. Although the 
details are somewhat different, the overall behaviour is the same 
as described above. In each case, r] A s wind penetrates deep into 
the acceleration zone of ? 7 B ’s wind and comes very close to the 
star at periastron. The key point is that the lower the value of M nA , 
the later the start of the cooling-transition phase and the earlier the 
recovery of the hottest X-ray emitting gas. In the Case B (Case C) 
simulation, the switch to the radiative-cooling regime occurs at 
c p R* 0.987 (0.991), with the hot 773 -shock apex fully recovering 
by 0 R 1.018 (1.014). The cause of this dependence on M, 1A is 
the decreased wind momentum of 77 a relative to that of 773 , which 
places the CD farther from r) B and allows its wind to obtain larger 
pre-shock velocities and lower post-shock densities. This increases 
the value of x since x oc u 4 ind and decreases the importance of 
radiative cooling in 773 ’s wind around periastron. 

The results of our 3D simulations are quite different from those 
of P 1 1 . The analysis in Appendix A2 can be used to help understand 
why. First, while the wind of 773 in our simulations is parametrized 
using a 7 S = 1 velocity law, Pll use standard CAK parameters for 
? 7 b that are more consistent with fi R 0.8. Our simulations also use 
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Figure 3. First three columns: Log density, temperature and total wind speed (left to right, respectively, cgs units) in the orbital plane at phases 0 = 0.98, 0.985, 
0.99, 1.00, 1.01 and 1.02 (rows) from the very small domain 3D SPH simulation of 77 Car assuming A/^ a = 8.5 x 10 ~ 4 Mq yr -1 (Case A). All plots show the 

inner ± la region. Axis tick marks correspond to an increment of 0.1 a « 1.54 au. Rightmost column: line plots of the total wind speed iyj ^ + Vy + v£) along 
the line of centres between 773 and from the simulation at each phase. Distance is measured from the surface of 773 in au. Note that all line plots extend only 
to r = 2 au from 773 in order to emphasize the WWC and pre-shock wind speeds. The stellar radius of 773 is R*,^ & 0.14 au. 
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Orbital Phase ( <p ) 


Figure 4. Cooling parameter / in fj B ’s wind as a function of orbital phase 
0 assuming the simulation parameters in Table 1 (Cases A— C) and PI 1 . 
The pre-shock wind speeds are computed using equations (A15)-(A19). 


a slightly larger R* for jy B (30 Rq versus 20 Rq). As illustrated 
in Fig. 6(a), even without RI. larger R , and /3 lead to substantially 
lower rj B wind speeds at distances < 1 .5 au from r] B . For our adopted 
R», a [5 = 1/2 would be required to roughly match the ?j B wind- 
velocity profile used by Pll. However, simply reducing R* would 
also result in a wind- velocity law that closely matches that of PI 1. 
Thus, the assumed R, of tj B can have as much of an influence on 
the overall wind-velocity profile as the /3 index. 


Inclusion of the velocity-altering effects of RI makes our pre- 
shock r/ B wind speeds even lower compared to those of P 1 1 . Setting 
1 1 = 0.8 in equation (A4) and substituting this into equation (A12), 
we can numerically integrate equation (A12) using the simulation 
parameters of Pi 1 . We can then compare the results to those for the 
parameters in Table 1 . To allow for a more direct comparison with 
Pll, we investigate here parameters for our Case B simulation with 
M Va = 4.8 x 10~ 4 Mq yr~' and L m = 3.3 x 10 5 Lq. 

Fig. 6(b) shows the results for several D corresponding to <j> = 0.5, 
0.96, 0.98 and 1 .0. At apastron, the difference between our expected 
total pre-shock ?j B wind speed and that using the parameters of Pi 1 
is relatively small, ~100kms~ 1 . This is true for most of the orbit. 
However, the difference between the two models becomes much 
more pronounced as periastron is approached. At <f> = 0.96, our pre- 
shock r/ B wind speed is ~440kms _1 lower than that of PI 1, while 
at </> = 0.98, it is ~815kms _1 lower. The largest, most important 
difference occurs at periastron, where for our model parameters )/ B 
is unable to drive a wind along the line of centres. r/ B can, however, 
drive a wind to a rather large pre-shock velocity of ~1920kms~ 1 
when using the parameters of PI 1 . 

PI 1 speculated that an increase in the wind velocities when mov- 
ing into periastron helped stabilize the WWC in their 3D simulations 
and prevented a collapse. While this may be true, we find that the R* 
and wind-velocity law assumed for ij B also have a big influence on 
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Figure 5. Same as the first three columns of Fig. 3, but for simulation Cases B (left three columns) and C (right three columns). 
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Figure 6. wind speed along the line of centres as a function of distance from ;j B . Left panel: standard /1-wind- velocity laws computed for the combinations 
of ft and R t that correspond to the ij B parameters of Table 1 (blue) and PI 1 (red). Included is the velocity profile for ft = 1/2 and R* = 30 Rq (dashed black 
line) that also roughly matches the velocity profile of used by PI 1 . Right panel: total q B wind speed i> 2 , tot via equations (A15)— (A19), using the parameters 
in Table 1 with M, jA = 4.8 x 1 0 Mq vr 1 . q = 15.22 (the Case B simulation, dashed lines) and the simulation parameters of PI 1 (solid lines), for binary 
separations D corresponding to (j> = 0.5 (black), 0.96 (blue), 0.98 (red) and 1.00 (green). Each curve is truncated at the radius of ram pressure balance />, given 
by equation (A19). There is no stable balance at </> = 1 using the Case B parameters, implying a collapse of the apex of the WWC region. 


determining whether there will be a cooling-transition phase in q B s 
post-shock wind around periastron. The larger R t and f) used in our 
simulations lead to a significant reduction in the pre-shock q B wind 
speed, which results in strong, rapid cooling of the post-shock q B 
gas around periastron. Using the parameters of PI 1, we compute a 
X < 1 in the post-shock q B wind only between rf> ^ 0.997 and 1 .003 
(Fig. 4), and a minimum value of x ^ 0.64 at periastron, implying 
that the reduction in q B ’s wind speed by q A is insufficient to cause 
)/ B 's wind to switch strongly to the radiative-cooling regime. This 
helps explain the absence of a WWC collapse in the 3D simulations 
ofPll. 

While radiative forces and inhibition are important factors that 
help determine whether a cooling-transition phase and WWC col- 
lapse occur in q Car at periastron, we agree with Pll that the key 
factor is radiative cooling. Earlier 3D SPH simulations of q Car 
that launched the stellar winds at their terminal speed and were 
isothermal (Okazaki et al. 2008) or adiabatic (M12; Madura & 
Groh 2012) did not show any cooling-transition phase or WWC col- 
lapse. Investigating further, we performed two very small domain 
( r = 1.5a) 3D SPH simulations with M nA — 8.5 x 10~ 4 Mq yr _1 , 
focusing on phases around periastron. One simulation includes 
ft = 1 velocity laws, RI effects and assumes adiabatic cooling, 
while the second simulation includes radiative cooling, but no ra- 
diative forces or inhibition (i.e. the winds are launched at terminal 
speed). 

We find that a cooling-transition phase only occurs in the radiative 
cooling simulation (left-hand panel of Fig. 7). However, because the 
winds are launched at terminal speed and there is no RI, the pre- 
shock q B wind speeds are much higher and the cooling transition 
occurs over a much narrower range of phases ( ij> ^ 0.998—1.002). 
In the adiabatic simulation, the hot apex of the q B shock survives 
throughout the entire simulation and there is no cooling transition 
phase (Fig. 7, right-hand panel). Due to the twisted WWC zone and 
embedding of q B in q A s dense wind at late phases {ij> sa 0.995), 
the most distant parts of the trailing arm of the WWC region still 
cool to T < 10 6 K, producing a temperature asymmetry between the 
leading and trailing arms in the outermost post-shock q B wind. This, 
however, is not a collapse or radiative cooling/inhibition effect, but a 
consequence of rapid orbital motion in a highly eccentric orbit. Our 
results thus suggest that radiative cooling of the post-shock q B wind 
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Figure 7. Slices showing log temperature in the orbital plane at periastron 
from very small domain SPH simulations of q Car that use radiative cooling 
with the stellar winds launched at terminal speed (left-hand panel) and 
adiabatic cooling with ft = 1 velocity-law winds and RI effects (right-hand 
panel). All plots show the inner ±0.25o region. 

is primarily responsible for any disappearance of hot X-ray emitting 
gas and/or WWC collapse in q Car during periastron passage. Ra- 
diative forces (e.g. ft = 1 wind-velocity laws and inhibition) assist 
the cooling transition by helping slow q B ’s wind prior to periastron, 
allowing the transition to start earlier and last longer compared to 
a situation in which the winds collide at near their terminal speeds, 
since x oc u 4 ind . 

3.1.2 The influence of radiative braking 

Our antigravity approach does not currently permit the implemen- 
tation of radiative braking in the SPH simulations (for details see 
Gayley et al. 1997; Russell 2013). P09 demonstrated that q A cannot 
brake q B ’s wind, irrespective of what CAK parameters they used. 
Kashi & Soker (2009) claim there should be a normal ram bal- 
ance between q Car’s colliding winds at periastron and no radiative 
braking. However, P09 found that q B can radiatively brake q A s 
wind if it comes extremely close to the stellar surface, e.g. within 
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~80Rq (~0.37 au) of q B . While under normal circumstances this 
means radiative braking would not occur before the two winds 
collide, it implies that radiative braking may be able to prevent 
the collision of dense q A wind material with the surface of q B at 
periastron. 

The two conditions required for radiative braking can be sum- 
marized as (1) d = d/d rb > 1 and (2) P = P l2 / P lb > d 2 , where 
d = D/R 2 is the binary separation scaled by the secondary’s stel- 
lar radius, P n = (M t Vi)/(M 2 v 2 ) is the primary/secondary wind 
momentum ratio, P rb = 4f)^d 2 b /(2 + fi) 2+ P for a velocity law with 
index f), and d lb is found by solving the transcendental expression 
d rb = 1 + (t/rb/f ) (1 -“>/( 1 +“) for standard CAK power index a and 
constant f given by 


(. L 2 /L l ) l/( ~ l - a '> 2 GMi 
(1 + 0!)a“A 1-o) R 2 y^o i 


(Gayley 1995). 


( 1 ) 


Condition (1) implies that a photospheric collision of star l’s wind 
with the surface of star 2 is prevented, while condition (2) implies 
that there is no ‘normal’ ram pressure balance. 

Using the values in Table 1 with M VA = 8.5 x 10~ 4 MQyr _1 , 
q — 22.28 and a = 2/3, we see that condition (1) is satisfied through- 
out q Car’s 5.54 yr orbit, which means that q B ’s radiation is capable 
of preventing q A 's wind from impacting q B , even at periastron. For 
most of the orbit, condition (2) is not satisfied and there is no ra- 
diative braking. However, for 0.99 < <j> < 1.01, Equations (A15) 
through (A19) indicate that there is no normal ram pressure balance 
along the line of centres because of the RI of q B ’s wind by q A . If 
there is indeed no ram balance between <j> = 0.99 and 1.01, then both 
(1) and (2) are satisfied and there will be radiative braking of q A s 
wind by q B . This is consistent with the findings of P09. The reason 
Kashi & Soker (2009) determined that there would be no radiative 
braking at periastron is because they neglected the effects of RI 
and orbital motion and assumed that P\ 2 is constant throughout the 
orbital cycle. As shown in Appendix A2, q Car’s wind momentum 
ratio is not constant and depends strongly on the binary separation. 
While our analysis shows that radiative braking may occur between 
(f) = 0.99 and 1.01 in q Car, this does not significantly affect our 
results or conclusions since the main effect of such braking would 
be the prevention of the collision of dense primary wind with q B ’s 
surface at periastron if/when the WWC zone ‘collapses'. Radiative 
braking, however, has important implications for models in which 
q B is assumed to accrete wind material from q A at periastron, such 
as that proposed by Kashi & Soker (2009), since the braking could 
likely prevent such accretion. 


3.2 Large-domain simulations 

Figs 8 and 9 show, respectively, the density and temperature in the 
orbital plane at five select phases from the r = 100a simulations. 
Results for slices in the xz and yz planes are in Appendix B2, while 
movies of the simulations are available in the online version of this 
paper (see Supporting Information). 

The dense arcs and shells of q A wind visible in the outer (>20a) 
regions of Fig. 8 indicate that the binary has already undergone 
multiple orbits. Narrow cavities carved by q B in q A s wind during 
each periastron passage also exist on the back (periastron) side of 
the system. The width of these cavities increases with decreasing 
M 1]A due to the change in the wind momentum balance. The larger 
cavities in Cases B and C also contain hotter q B wind material and 
extend farther in the +y direction (Fig. 9). 

Bordering the wind cavities on the periastron side of the system 
are compressed, high-density shells of q A wind that form as a re- 


sult of the rapid WWC. These shells flow outward at v > u 00> ^ A , 
accelerated by the collision with q B ’s higher velocity wind. The 
speed of the shells increases with decreasing M, A , while the den- 
sity and thickness of the shells decrease with M nA . The innermost 
shell closest to the stars has the highest outward velocity due to the 
more recent collision with q B ’s wind. The outer shells move slower, 
at approximately v XtJ1A , due to the gradual increase in the amount 
of swept-up ?jA-wind mass in each shell. 

The counter-clockwise direction of orbital motion in the xy plane 
also produces a noticeable difference in the thickness of the shells 
and width of the q B wind cavities located in the +y and — y directions 
(Fig. 8). The shells of q A wind that lie in the — y direction are 
thinner and more compressed, while the q B wind cavities are wider. 
This size difference becomes more apparent as the value of M rlA is 
lowered, with Case C having the thinnest shells and widest cavities 
in the — y direction. 

Following periastron, q B returns to the apastron side of the sys- 
tem, its wind colliding with the q A wind that flows in the +x, —y 
direction. At </> rs 3.1, the arms of the WWC region become ex- 
tremely distorted by orbital motion and the leading arm collides 
with the old trailing arm from before periastron, forming another 
dense, compressed shell of q A wind that flows in the +x and — y 
directions (rows d and e of Fig. 8). The overall stability of this 
expanding shell depends on the shock thickness (Vishniac 1994; 
Wiinsch et al. 2010; Pll). Portions of the shell moving in the — y 
direction appear to be the most stable due to the large amount of q A 
wind that borders it in this direction. The shells are also more stable 
and remain intact longer the higher the value of (i/, A . In Case A, the 
upper part of the shell expanding in the +x, +y direction does not 
start to fragment until <^> fa 3.34. In Case C, the shell is completely 
disrupted in the +x, +y direction by </> 3. 1, the wind of q B having 

ploughed through it, causing it to separate from the leading arm of 
the WWC zone (see rows d and e of Fig. 8, and row e of Figs 1 and 
Bl). This separation produces a pair of dense ‘arcs’ of q A wind on 
the apastron side of the system. Multiple pairs of arcs, formed dur- 
ing previous orbital cycles, are visible in Fig. 8. They are also quite 
spatially extended; by the time the system is back at periastron the 
arcs formed during the previous periastron are up to 80a 1235 au 

from the central stars. 

As the dense arcs flow outward, the wind of q B collides with 
and propels them at a speed slightly greater than v 00iTlA into the low- 
density cavity created during the previous orbital cycle (Fig. 8). This 
produces another shock located ~20 a from the central stars (along 
the +x axis at <fi = 0.5) that heats the outer q B wind to T > 10 7 K 
(Fig. 9). The arcs also expand as they move outward. Since they 
are bordered on both sides by low-density, high-speed q B wind, the 
arcs eventually become unstable and fragment, mixing with the sur- 
rounding q B wind material. The lower the value of M, 1A , the thinner 
and less dense the expanding arcs, and the sooner they fragment 
and mix with the q B wind. An example of this fragmentation and 
mixing can be seen in the outermost arc/shell located at the bottom 
of the panels for Case C in Fig. 8. 


4 IMPLICATIONS FOR VARIOUS 
OBSERVATIONAL DIAGNOSTICS 

The SPH simulations show that for the range of explored values, a 
factor of 2 or more decrease in M r]A alters significantly the time- 
dependent 3D density, temperature and velocity structure of q Car’s 
spatially extended primary wind and WWC zones. Therefore, any 
such drop in M r]A should result in detectable differences in numerous 
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Figure 8 . Slices in the orbital plane from the large-domain 3D SPH simulations of 77 Car (Table 1) at orbital phases 0 = 2.5 (apastron), 2.7, 3.0 (periastron), 
3.15 and 3.30 (rows, top to bottom). Columns correspond to the three assumed 77 a mass-loss rates. Colour shows log density in cgs units. The spherical 
computational domain size is r = 100a « 1545 au ~0.67 arcsec ( D = 2.3 kpc). Axis tick marks correspond to an increment of 10 a & 155 au ~0.067 arcsec. 
The orbital motion of the stars is counter-clockwise. 77 a is to the left and tjb is to the right at apastron. 
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Figure 9. Same as Fig. 8, but with colour showing log temperature. 


observational diagnostics (compared to earlier cycles). Guided by 
the SPH results, we qualitatively discuss below some of the expected 
changes and suggest observations to perform in order to further test 
whether is decreasing. We focus here on changes that are 


likely to be most important. Detailed quantitative modelling of the 
vast array of multiwavelength observations of 17 Car using the 3D 
simulations is beyond the scope of this paper and is left for future 
work. 
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4.1 Effects on the amount of material in line of sight 

Ml 2 tightly constrained, for the first time, the 3D orientation of 
)/ Car’s binary orbit using a 3D dynamical model for the broad, 
spatially extended [Fe hi] emission observed by the HST/ST1S (Gull 
et al. 2009). M12 find that the observer’s line of sight (LOS) has an 
argument of periapsis cu tv 240° to 285°, with the binary orbital axis 
closely aligned with the Homunculus polar axis at an inclination 
i sa 130° to 145° and position angle on the sky PA. sa 302° to 327°, 
implying that apastron is on the observer’s side of the system and 
that tj b orbits clockwise on the sky (see rightmost panel of Fig. 1 0). 
This orientation and clockwise motion of rj b are consistent with all 
known observations of rj Car to date (see the discussion in M12 and 
Madura & Groh 2012). 

Assuming the orbital orientation in Fig. 10, a factor of 2 or more 
drop in M nA , for the range of explored values, would change signif- 
icantly the time- variable amount of material in LOS to the observer. 
For example, if M, 1A = 8.5 x 1CL 4 Mg yr , then the half-opening 
angle of the WWC zone (^55°) is quite close in value to the inclina- 
tion angle of the binary away from the observer (^48°). Therefore, 
the observer’s LOS lies very near the edge of the WWC cavity 
(see bottom-left panel of Fig. 10), with portions of the dense, outer 
interacting winds passing through our LOS. The column density 
is thus expected to be greater for sightlines to the observer than 
to most other parts of the system, especially towards the Weigelt 
blobs, which lie in/near the orbital plane (Davidson et al. 1997; 
Smith et al. 2004a) surrounded by the low-density wind from fjB- 


As discussed in G12a, this may resolve the long standing puzzle 
of why our direct view of the central star suffers more extinction 
than the Weigelt blobs located only 0.3 arcsec away (Davidson et al. 
1995; H01). Furthermore, as suggested by Smith (2010) and G12a, 
dust formation in the densest parts of the outer interacting winds 
might explain the peculiar ‘coronograph’ described by H06. 

A decrease of less than a factor of 2 in M nA is enough to increase 
the opening angle of the WWC cavity and cause the observer's 
LOS to no longer intersect portions of the WWC region (bottom 
row of Fig. 10). For A = 8.5x 10~ 4 Mq yr -1 , the difference in 
angle between the observer’s LOS and the dense wall of the WWC 
zone is ~5°-10°, whereas for 4.8 x 10~ 4 Mq yr~ ! , the difference 
is ~15°-20°. Fig. 10 shows that the lower the value of the 
more the observer’s LOS samples only low-density wind. The 
column density in LOS to most other directions on the apastron 
side of the system, e.g. towards the Weigelt blobs, is not nearly as 
strongly affected by decreases in M nA . Interestingly, since 1998, 
the Weigelt blobs have not brightened significantly, even though the 
central source has brightened by one stellar magnitude in the optical 
(Davidson 2005; Gull et al. 2009). A change in rj caused by a slight 
decrease in M nA might explain this since the larger opening angle 
would decrease the amount of extinction in LOS to the observer, 
causing an apparent increase in the brightness of the central source 
while leaving the Weigelt blobs unaffected due to their location in 
the orbital plane. Additional observations are needed in order to 
determine that such a change in opening angle occurred and is due 



Figure 10. Geometry and location of the colliding stellar winds relative to the observer and the central stars, and incorporating the apparent position of two 
Weigelt blobs, at an orbital phase near apastron. Top row of three panels: slices showing log density in the orbital plane from the 3D SPH simulations for 
the three different values of M ljA . The locations of the stars, wind cavity and Weigelt blobs are indicated by arrows. The dashed line and eye indicate the 
projected line of sight to the observer. Axis tick marks correspond to an increment of 10a 155 au ~0.067 arcsec. Bottom row of three panels: same as top 

row. but for slices in the xz plane. Far right panel: illustration of rj Car’s binary orbit (upper left inset, yellow) on the sky relative to the Homunculus nebula 
assuming the binary orientation derived by M12 ( i = 138°, or = 263°, and PA. = 317°). The +z orbital axis (blue) is closely aligned with the Homunculus 
polar axis in 3D. r/B orbits clockwise on the sky relative to ij ,\ (black arrows in inset) and apastron is on the observer’s side of the system. The semimajor (+x) 
and semiminor (+y) axes are indicated by the red and green arrows, respectively, and North is up. Included for reference are insets showing an HST WFPC2 
image of the Homunculus (top right, Credit: NASA, ESA and the Hubble SM4 ERO Team) and a 1 arcsec x 1 arcsec HST / ACS F550M image (lower right) 
showing the location of the central stars (A) and the bright Weigelt blobs (B, C, D). 
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solely to a decrease in M, A and not, say, an increase in the M and/or 
Vcc of riB- 

4.2 Effects on X-rays 

4.2. 1 Column density and the X-ray minimum 

Changes to the wind momentum ratio will change the WWC open- 
ing angle and alter the distribution of the shock-heated gas, which 
should have observable consequences in the X-ray band. An in- 
crease in the WWC opening angle of ~10°-20° should cause a 
delay of ~20-40d in the observed sharp rise in column density 
(Ah, Hamaguchi et al. 2007) that occurs prior to periastron. The 
decrease in Ah seen after periastron should also begin ~ 20-40 d 
earlier. When combined with the reduced duration of the cooling- 
transition phase of jjb's post-shock wind, the net result should be a 
shorter X-ray minimum during periastron passage for lower values 
of M, 1a . This simple analysis reveals that for the range of explored 
M, a , the total X-ray minimum should be ~40-80d shorter than 
that of a previous orbital cycle. 1 * 

One mystery regarding rj Car’s X-ray light curve is why no signif- 
icant change was observed in 2009 just prior to the X-ray minimum. 
Based on the simulations, the start of the 2009 minimum should have 
also been delayed if M, 1A dropped by a factor of 2 or more. Yet, only 
the timing of the X-ray recovery changed. One might argue that this 
is evidence of a sudden significant (e.g. > factor of 2) change in 
M nA during the 2009 event, with the apparent observed decrease in 
X-ray hardness between 2004 and 2009 due to a smaller, gradual 
decrease in M nA that was not enough to change the timing of the 
start of the 2009 X-ray minimum. Since the recovery from the 2009 
event, the 3-10 keV X-ray flux has been ~20 percent below the 
level seen in the previous two cycles (Corcoran 2013). A decrease 
in flux outside of the minimum is expected if M J]A decreased. How- 
ever, a drop in M, 1A by a factor of 2 should result in a flux decrease 
of about a factor of 2 (Russell 2013). The implication is that, if 
M 1]a decreased around the 2009 event, either (1) the change was by 
a factor < 2 or (2) the decrease occurred suddenly and by a fac- 
tor > 2, but immediately after the 2009 event, M IjA quickly returned 
to its ‘normal’ value, resulting in an overall X-ray flux decrease of 
~20 percent at times away from periastron. 

Furthermore, the 2003 X-ray minimum did not start earlier, which 
seems to rule out any large change in M, 1A prior to mid-2003. The 
first significant sign of a possible recent change in Car was the 
early recovery of the 2009 X-ray minimum. Up until that point, 
most differences in ij Car were viewed as typical cycle-to-cycle 
variations and/or caused by the continuous decrease in the amount 
of circumstellar extinction in LOS (Martin, Davidson & Koppelman 
2006a; Smith 2010). Detailed analysis by CIO of r] Car’s X-rays 
over the last three cycles revealed later that the X-ray hardness 
and flux had decreased slightly between 2004 and 2009. The X-ray 
spectrum was also softer 100 d after the 2009 minimum than it was 
100 d after the 2003.5 minimum, with the E > 3keV flux lower as 
well. Moreover, it now appears that a high M, A of ~10~ 3 Mg yr~ J 
is necessary in 3D hydrodynamical simulations to reproduce the 
extended RXTE X-ray minima in 1998 and 2003.5, and the X-ray 
flux level observed between events (Russell 2013; Russell et al., 
in preparation). Therefore, the X-ray data suggest that if there was 
a recent large decrease in M, 1A , it occurred sometime around/after 

1 For detailed quantitative modelling of r] Car’s X-ray observations that 

support this result, see Russell (2013). 


2004. The lack of any change in the Hen 7.4686 emission before 
the 2009 event supports this (see Section 4.3). 

If the extended X-ray minima of 1998 and 2003.5 were the result 
of a long cooling-transition phase, perhaps the shorter minimum in 
2009 was the result of some change in the effectiveness of radiative 
cooling in tjb’s post-shock wind. A modest decrease in M Va of 
~20 percent (or a bit more) might have moved the WWC apex far 
enough from ijb to decrease the effectiveness of radiative cooling, 
but a factor of 2-4 drop does not appear to be fully consistent 
with the current X-ray data. On the other hand, if M, 1A is indeed 
decreasing and continues to do so, the 2014 X-ray minimum should 
be even shorter than that of 2009. Continued decrease of M Va would 
eventually lead to an even shorter X-ray minimum. 

4.2.2 The CCE component 

During the deep minimum, the X-ray emission is dominated by the 
‘Central Constant Emission’ (CCE) component (Hamaguchi et al. 
2007, hereafter H07). Chandra data indicate the CCE is pointlike 
(< 2000 au for D = 2.3 kpc). The X-ray luminosity of and column 
density to the CCE component suggest that it is a stable (over the 
course of the deep minimum), luminous X-ray source that exists 
beyond the primary wind in front of the WWC apex at periastron 
(H07). The large-domain SPH simulations (Figs 8, 9 and B5— B8) 
show that, at periastron, the collision of ^ B ’s wind with the arcs 
of compressed a wind material located on the apastron side of 
the system could produce significant amounts of hot (>10 7 K) gas 
in front of the primary at distances <2000 au in size. The simula- 
tions further show that this emission varies slowly on time-scales of 
months, and so this faint, hot gas is a logical candidate for the source 
of the CCE. Any decrease of M VA will propagate into the outer wind 
structures of r) Car, affecting the thermalization of the CCE com- 
ponent and/or the absorption to the CCE in LOS. Measurements of 
the X-ray emission measure, column density in LOS and ionization 
state of the CCE component during Car’s next periastron in 2014 
would thus be useful. When compared to the same data from the 
2009 event, the 2014 observations could help confirm or rule out a 
recent factor of 2-4 drop in M, 1A . 

4.2.3 X-ray emission line profiles 

X-ray emission line profiles provide important information on the 
temperature and density distribution of the hot, post-shock ijb gas 
flowing along the boundary of the WWC zone (Henley et al. 2008). 
They can also be used to constrain the opening angle of the WWC 
cavity and the orientation of the system relative to our LOS (Henley 
et al. 2008, PI 1 ). Chandra high-energy grating spectra show that 
strong, narrow emission lines are relatively unshifted during most 
of r) Car’s orbit. lust prior to the X-ray minimum (cp & 0.98-0.99), 
the strongest lines become increasingly blueshifted and broader, 
reaching velocities 2000 km s~* as the flow of shocked ;; B wind 
along the trailing arm of the WWC zone swings rapidly past the 
observer’s LOS (Henley et al. 2008; Corcoran 2011). 

The small-domain simulations show that an increase in the WWC 
opening angle causes the hot trailing arm of the WWC zone to cross 
the observer’s LOS (with a> & 260°) at later orbital phases as the 
system approaches periastron. The observer will thus not see hot, 
high-velocity post-shock tj b gas until much closer to periastron. 
This should result in a detectable phase delay in when the X-ray 
lines become blueshifted and broader. The magnitude of this delay 
will be proportional to the increase in WWC opening angle. 
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4.2.4 X-ray flaring 

ri Car’s RXTE light curve also exhibits short-term rapid varia- 
tions in brightness or ‘flares’ (Corcoran et al. 1997; Moffat & 
Corcoran 2009; CIO). These flares become more frequent and 
shorter in duration as periastron is approached. Moffat & Corcoran 
(2009) presented a series of models for explaining the flares, their 
preferred model consisting of large stochastic wind clumps from 
riA that enter and compress the hard X-ray-emitting WWC zone. 
As an alternative, they considered that the flares are the result of 
instabilities intrinsic to the WWC zone. Pll indicate that such in- 
stabilities, especially if combined with a clumpy LBV wind, may 
be a viable explanation for the X-ray flares. The SPH simulations 
suggest that a decrease in the value of M, 1A leads to a more un- 
stable, clumpy WWC region. A decrease in M r]A would therefore 
result in a higher flare frequency, especially as periastron is ap- 
proached. The relative duration and amplitude of the flares may also 
change, although detailed numerical modelling is required to test 
this. 

4.2.5 Changes to tie ’s wind parameters? 

An increase in the opening angle of the WWC zone combined with a 
shortened cooling-transition phase during periastron passage might 
explain the shortened duration of the 2009 X-ray minimum. How- 
ever, these effects may not necessarily be a result of a drop in M, A . 
A comparable increase in the M and/or of )j B ’s wind would also 
increase the WWC opening angle. One argument against a recent 
increase in the i>oo of tie’s wind is that such a change should increase 
the X-ray hardness at all phases. However, the X-ray hardness ap- 
pears to have declined somewhat between 2004 and 2009 (CIO). A 
significant increase in the M of rj B is harder to quantify since the 
densities in the post-shock ?j B wind would increase, affecting the 
X-ray emission measure. Recent flux measures are lower compared 
to previous cycles, suggesting a decrease in the emission measure 
of the shocked ;; B wind. One argument against an increase in ?? B ’s 
M is that such an increase should strengthen the effectiveness of ra- 
diative cooling in the post-shock wind at phases close to periastron. 
One would then expect a longer cooling-transition phase that leads 
to a longer minimum, as opposed to the observed shorter minimum. 
Finally, a decrease in stellar radius or change in the velocity-law 
profile (e.g. ft) of t]e could also allow its wind to reach higher pre- 
shock velocities without drastically altering Voo, moving the WWC 
zone closer to ija and again decreasing the effectiveness of radiative 
cooling during periastron passage. Such changes could in principle 
produce a shortened X-ray minimum compared to previous cycles, 
although the physical reasons for a sudden alteration of rj B ’s stellar 
parameters seem inexplicable. Numerical simulations investigating 
tie’s stellar/wind parameters are needed in order to determine if 
such changes can explain the X-ray observations. 

4.3 The He II 74686 emission 

Somewhat related to the X-ray emission is the peculiar behaviour 
of the He ii 74686 emission around periastron (Steiner & Damineli 
2004; Martin et al. 2006b; Mehner et al. 2011; Teodoro et al. 2012). 
The He n 74686 line intensity increases suddenly by a factor of 10 
just before phase 1.0 and then drops sharply to zero, after which 
it recovers to a second peak in intensity before declining back to 
zero again. A highly luminous source of He + -ionizing photons with 
hv > 54.4 eV is required for the formation of this line. The necessary 
soft X-rays are thought to come from the shocked t/a wind, while 


the He n 74686 emission is hypothesized to originate in the 77 A wind 
just before the WWC and/or in the cold, dense post-shock tia wind 
(Martin et al. 2006b; Abraham & Falceta-Gon§alves 2007; Mehner 
et al. 2011; Teodoro et al. 2012). Earlier alternative explanations 
in which the He ii 74686 arises in the wind of t/e can be found in 
Steiner & Damineli (2004) and Soker & Behar (2006). 

Mehner et al. (2011) argue that the Hen emission behaved dif- 
ferently during the 2009 spectroscopic event and included a second 
intensity peak just after periastron that was not present during the 
2003 event or delayed. However, Teodoro et al. (2012) claim that the 
two events were not drastically different. In any case, both Mehner 
et al. (2011) and Teodoro et al. (2012) support a scenario in which 
there is a WWC ‘collapse’ during periastron that helps explain the 
behaviour of the He ii 74686 line. A WWC disruption seems neces- 
sary to explain the fact that the He n behaviour is remarkably similar 
for both reflected spectra from the polar region of i] A and a direct 
view of the central stars (Stahl et al. 2005; Mehner et al. 201 1). 


4.3.1 Formation of the Hen 74686 emission in t) Car 

Before elaborating on how the He ii emission should change if there 
is a large drop in /i/, A , it is important to clarify how/where the He u 
emission forms in i] Car, and why the equivalent width (EW) be- 
haves the way it does around periastron (see e.g. fig. 4 of Mehner 
et al. 2011 or fig. 3 of Teodoro et al. 2012). Some of the ideas below 
have been discussed in Martin et al. (2006b); Mehner et al. (2011) 
and Teodoro et al. (2012), to which we refer the reader for details. 
However, we present several new concepts that are crucial for un- 
derstanding the formation of the observed He n 74686 emission and 
its dependence on M^ A . 

As discussed by Martin et al. (2006b); Mehner et al. (201 1) and 
Teodoro et al. (2012), the low-density pre-shock wind of t/e, the 
acceleration region of the wind of tie and the post-shock wind of ;j B 
can all easily be excluded as the source of the He ii 74686 emission, 
based simply on energy and density requirements. While it is often 
assumed that the Hen emission arises in the dense wind of ?; A , 
two key remaining questions are (1) why is strong Hen emission 
only observed between tj> ~ 0.98 to 1.03 and (2) how/where do the 
required He + ions form? 

According to fig. 9 of H01, for r > 10R* iI)a , helium is neutral in 
the dense wind of ?; A (assuming M rlA = 10~ 3 Mq yr _1 ). Therefore, 
under normal circumstances, there can be no Hen emission for 
most of ti Car’s cycle since there is no He + near the WWC zone 
to be ionized. Moreover, for most of the orbit, soft X-rays formed 
in the WWC zone cannot penetrate the dense primary wind that 
lies between ;; B and the small He + zone at r < 10 R*, nA . There is 
no other known source of photons capable of ionizing the inner ri A 
He+ zone. One might try to argue that radiation from the WWC 
shocks ionizes the neutral helium in t)a ’ s wind located near the 
WWC region. However, the WWC shocks are unsuitable for this 
since the 3000 km s~' shocked wind of ?7 B produces mainly hard 
X-rays with energies above 500 eV, while the 420kms~ 1 wind of 
t] A produces photons with energies of the order of a few hundred 
eV. Both are highly inefficient for ionizing neutral helium, and only 
the shocked wind of t/a is suitable for providing photons to ionize 
He+. 

The only other possible source of neutral-helium ionizing photons 
at times away from periastron is t/e ■ As noted by Nielsen et al. (2007) 
and Humphreys, Davidson & Koppelman (2008), during most of the 
binary orbit, ?? B may photoionize neutral helium in a small region of 
the pre-shock tia wind located just beyond the WWC zone (similar 
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to fig. 12 of Humphreys et al. 2008)3 Neutral helium in portions 
of the colder, denser post-shock i;a wind would also be ionized by 
tj B . The He + in the dense post-shock q A wind will very quickly 
recombine, producing He i emission. Some of the He + in the lower- 
density, pre-shock q a wind would also recombine and create He l 
emission (Humphreys et al. 2008). However, we suggest that ifq B 
can photoionize neutral helium in the pre-shock wind of q A , some 
of the resulting He + will be additionally ionized by soft X-rays 
created in the colliding q A wind, producing a very small amount of 
He n 24686 emission at times away from periastron. The amount of 
He n emission is likely to be small because of the lower densities in 
the pre-shock q A wind near the WWC region, and because of the 
likely very small volume of He + created by q b’s radiation. This is 
consistent with the observations and explains why little to no He ii 
2.4686 emission is detected during most of q Car’s 5.54 yr cycle 
(Steiner & Damineli 2004; Stahl et al. 2005; Martin et al. 2006b; 
Mehner et al. 201 1 ; Teodoro et al. 2012). 

Earlier proposed scenarios for the He n 2.4686 emission also have 
difficulty explaining the emission peaks observed just prior to and 
after periastron. Here, we present a new, simple explanation for 
these peaks, fig. 9 of H01 shows that He + exists in ?; A ’s wind only 
between ~2.4R tJ1A and 10/?*.,, A (0.7-3 au). The densities in this 
region (~10 u -10 13 cm~ 2 3 ) are capable of explaining the observed 
He ii emission if the He+ can be ionized by soft X-rays. If we assume 
that the 54-500 eV photons needed to ionize the He + come from 
the post-shock q A wind, we find that for most of q Car’s orbit, the 
WWC zone is too far from q A s He + zone to produce any Hen 
emission (Fig. 1 1). 

Plotting the separation between q A and the WWC apex as a 
function of phase assuming the parameters in Table 1, we see that 
only for 0.986 < <f> < 1.014 does the WWC apex penetrate into 
the He+ zone of q A (bottom panel of Fig. 11). Therefore, the soft 
X-rays from the WWC region can ionize the He + only during these 
times. This is almost exactly the same range of spectroscopic phases 
that strong Hen 24686 emission is detected. Furthermore, this is 
roughly the same range of phases that the post-shock q B wind is in 
the radiative-cooling regime. During this time, the post-shock q B 
will also produce photons capable of ionizing He + . 

Martin et al. (2006b) and Mehner et al. (2011) speculated that 
high densities in the post-shock q B gas around periastron could lead 
to a switchover to the radiative-cooling regime, which would make 
the WWC zone unstable and collapse, and lead to the increased 
generation of the soft X-rays required to form He n. According to 
our SPH simulations, a complete WWC collapse is not required 
for the additional generation of soft X-rays to ionize He + . A mass- 
ejection event, accretion by q B and/or changes to a latitudinally 
dependent primary wind (Martin et al. 2006b; Soker & Behar 2006; 
Mehner et al. 2011) are also not required. Instead, we suggest that 
the strong Hen 24686 emission arises in the inner He + zone of 
q A & dense wind and that the first peak in EW observed just before 
the spectroscopic event (0 m 0.98-0.995) is due to a combina- 
tion of two important effects; (1) penetration of the WWC region 
into the He + zone of q A that allows soft X-rays generated in the 
WWC zone to ionize the He + and (2) a switchover of the post- 
shock wind of q B to the radiative-cooling regime that results in 
additional soft X-rays that can ionize the He + and produce extra 


2 We say ‘may’ because it is possible that the dense post-shock q A wind will 
prevent q B from ionizing neutral helium in the pre-shock q A wind, in which 
case there would be no He n 24686 emission. Detailed 3D radiative transfer 

simulations are needed to test this though. 




Orbital Phase (0) 


Figure 11. Plot of the separation between r) A and the apex of the WWC 
zone as a function of orbital phase 0, assuming the parameters in Table 1 
with M, a = 8.5 x I (r 4 Mq yr -1 , q = 22.28, and using Equations (A15) 
and (A19) to determine the location of ram pressure balance ('/>, ). Apastron 
occurs at 0 = 0.5 and periastron at 0 = 0, 1. The solid-red horizontal line 
indicates the radius in the wind of q A where, based on fig. 9 of H01, helium 
transitions from being neutral to singly ionized. The dashed-red line marks 
the radius in q A ’s wind where helium transitions from being singly to doubly 
ionized (r ~ 0.7 au). The lower plot is a zoom of the dashed-box region in 
the upper plot (0.95 < 0 < 1.0). 


Hen emission. The fact that the timing of both matches almost 
perfectly with the observed first He n peak provides strong support 
for our scenario. The sharp drop in Hen EW at <j> ^ 1.0 may be 
due, at least in part, to a physical collapse of the WWC zone at 
periastron. 

We further suggest that the second, smaller emission peak in 
He ii 24686 that occurs ~+20 d after 0 = 1 .0 represents the period 
just after perisatron when q B returns to our side of the system and 
begins to carve a new cavity within ?7 a’s wind. During this time, 
the WWC zone may still be highly unstable, but two important 
facts will remain. Namely, the WWC region will be within the He + 
zone of q A and the post-shock wind of q B will be in the radiative- 
cooling regime and supplying additional soft X-rays to produce 
He ii. The decreased height of the second He n peak could be due 
to a combination of instability in the WWC zone and increased 
optical depth in LOS caused by the dense q A wind that flows in the 
apastron direction during periastron passage. The He ii behaviour 
should look relatively similar for both a direct view of the central 
source and reflected spectra corresponding to high stellar latitudes, 
with some possible changes in the line profiles. 
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4.3.2 Behavior of the Hen A4686 with decreasing M VA 

If M nA is lowered to ~2.5 x 1 0 -4 Mq yr - 1 , keeping the other stellar 
and wind parameters constant, ID cmfgen models show that r) A ’s 
He + zone increases greatly and extends to r ~ 120 au. As a result, 
the WWC zone is within the He + region for the entire 5.54 yr cycle. 
Soft X-rays generated in the post-shock r/ A wind should ionize the 
He + and produce significant, detectable He ii emission at all phases, 
with the exception of a possible short, sharp drop in He n EW at 
periastron. This emission, though, is not observed. 

If is lowered to ~5.0 x 10 -4 Mq yr -1 while keeping the 
other parameters constant, cmfgen models show that r) A s He+ zone 
extends to r ~ 7.5 au, more than a factor of 2 larger than in the 
1 0 4 Mq yr -1 model. Combined with the movement of the WWC 
apex closer to tj A , this should lead to a significant change in the 
timing of the He ii emission peaks. The WWC apex is predicted to 
penetrate the He + zone at 0 ~ 0.94, or ~85 d earlier. The first He n 
peak should therefore also appear ~85 d earlier. While this first 
peak may appear earlier, its initial amplitude may not be as large. 
This is because a lower M nA produces a shorter cooling-transition 
phase in jj b ’ s post-shock wind, resulting in the generation of less 
soft X-rays to ionize the He + . Therefore, an initial, broad He ii peak 
could begin at f ~ 0.94, created when the WWC apex penetrates 
the He + zone and the soft X-rays formed in the post-shock >i A wind 
produce He n emission. Another sharper peak could then occur near 
< p ~ 0.99 when tjb’s wind switches to the radiative-cooling regime 
and produces an additional burst of soft X-rays to further ionize 
the He + and create additional He n emission. Since the WWC zone 
will be more stable in a lower M, 1A situation, the sharp drop in He n 
EW near cf> = 1.0 may be delayed, while the recovery to a second 
He ii peak after periastron may be sped up. Because the WWC zone 
will remain in the He + zone longer after periastron, the height and 
duration of the second He ii peak should both increase. The time 
taken for the Hen to decline back to near zero EW should also 
increase. 

Based on these results, it is clear that for the explored parameter 
space, any change in M VA of a factor >2 should lead to an equally 
significant change in the behaviour of the He u A4686 EW. Inter- 
estingly, as shown by both Mehner et al. (2011) and Teodoro et al. 
(2012), the first Hen peak during the 2009 event occurred almost 
exactly 2023 d after the 2003.5 event. This and the fact that the sharp 
drop in X-rays also occurred nearly 2023 d after the 2003 periaston 
very strongly suggest that there was no large decrease of M nA before 
2009. As the available observations make it unclear whether a sec- 
ond He ii peak occurred during the 2003 event, there is a possibility 
M rlA dropped during the 2009 event. However, such a drop would 
likely have to be by a factor of < 2 (in contrast to CIO), otherwise, 
there would be a detectable increase in the amount of He ii L4686 
emission at phases far from periastron. Such emission is currently 
not seen. If M r]A has decreased, or continues to do so, changes in 
the behaviour of He ii X46S6 during the next spectroscopic event in 
2014 should be clearly visible. 


4.4 Effect on line profiles of broad wind-emission features 

4.4.1 Expected changes 

Using 2D radiative transfer models, G12a showed that the WWC 
cavity significantly affects the formation of Ha, H ft and Fe n lines in 
A ’s dense extended wind, even at apastron, by reducing the amount 
of P-Cygni absorption in LOS. Ionization of portions of rj A ’s ex- 
tended wind by tj b also affects the shape and strength of these stellar 


wind lines (G12a; Groh et al. 2012b, hereafter G12b). The WWC 
cavity furthermore induces a latitudinal and azimuthal dependence 
in these line profiles at orbital phases around apastron, explaining 
the stronger P-Cygni absorption observed in spectra scattered off 
of the poles of the Homunculus (G12a,b), which was originally 
interpreted as being caused by rapid stellar rotation (Smith et al. 
2003; Mehner et al. 2012). Observed latitudinal changes in the line 
profiles during periastron, specifically the appearance of strong P- 
Cygni absorption in LOS, is explained by the flow of dense primary 
wind towards the observer for a brief period when r) B is behind r] A 
at periastron (G12a,b; Madura & Groh 2012, see also rows c and d 
of Figs 2 and B 1). 

Since the broad wind-emission features of Ha, H/l and Fe u are 
so greatly affected by the WWC cavity, a change to M, lk that al- 
ters the WWC opening angle should have significant observable 
consequences for these and similar lines. We can make some pre- 
dictions using the SPH simulations and the previous spectral models 
of H01, H06 and G12a,b. Fig. 11 of H06 shows that, keeping the 
other parameters (e.g. R„, L„, etc.) fixed, as M, A is lowered, the 
Balmer P-Cygni emission and absorption in LOS weaken, with 
the absorption eventually disappearing. The Fen emission lines 
weaken as well, but emission lines of He i strengthen. Eventually, 
for low enough values of M^ a (~2.5 x 10 -4 Mq yr -1 ), lines of Niii 
and He n X4686 come into emission (H06). 

It is important to note that the results of H01, H06 were de- 
rived from ID spherically symmetric radiative transfer models and 
do not take into account the presence of the WWC cavity or r] B . 
Inclusion of the WWC cavity will further reduce the amount of 
Balmer P-Cygni absorption in LOS (G12a). The effect of the cav- 
ity on the absorption for a direct view of the central source may 
be difficult or impossible to detect though since for low values of 
M, 1a (~2.5 x 1CL 4 Mq yr -1 ), the P-Cygni absorption is sometimes 
completely absent in the ID models, depending on the assumed 7?*. 
There is no P-Cygni absorption in Ha in the models only if R, JJA is 
small (~60Rq). Assuming L„ = 5 x 10 6 Lq, a model with 60 Rq 
has a temperature at the hydrostatic core of T* ^ 35 000 K, while 
the model with 300 Rq has & 15 000 K. This results in a big 
difference in the ionizing fluxes between models, which ultimately 
determines the amount of Ha absorption. 

However, the increase in cavity opening angle will alter the lati- 
tudinal and azimuthal variation of the P-Cygni absorption (G12a). 
The photoionization cavity created by j/ b will also increase in size. 
Based on Fig. B 1, a factor of 2 or more drop in M ItA should reduce 
the r) A wind density, and thus the amount of P-Cygni emission and 
absorption in Ha, over the stellar poles (see also H06). A factor 
of 4 drop in M 1]A might result in almost no P-Cygni absorption in 
Ha over the poles of r) A at phases far from periastron (assuming 
that the dense, post-shock ;; A wind is ionized). The models predict 
that any larger drop in M 11A will increase the opening angle of the 
WWC cavity in such a way that it starts to bend towards t] A , likely 
resulting in little to no P-Cygni absorption in Ha at high stellar 
latitudes during most of Car’s cycle. 

The SPH simulations also suggest that a drop in M, A should re- 
sult in the detection of less P-Cygni absorption in LOS in Ha and 
similar lines during periastron passage. This is because, as shown in 
Figs 1 and B 1, as M,, A is decreased, the volume of t] A wind that flows 
towards the observer during periastron also decreases, producing a 
thinner, less-dense wall of rj A wind between the star and the observer 
after periastron. If M I]A is low enough (~2.5 x 10 -4 Mq yr -1 ), for 
R, = 60 Rq, there may be no increase in P-Cygni absorption in 
LOS during periastron since >j A will be able to fully ionize hydro- 
gen in this outflowing primary- wind region (H06). The decreased 
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Figure 12. Comparison of the spectrum of ij Car observed with //.ST/ST IS in 2001 Apr 17 (ip = 10.603; black solid line) with ID cmfgen models assuming 
Voo = 420 km s 1 with M rjA = 8.5 x 1 0 4 Mg yr - 1 (blue dashed line) and 2.5 x 1 0 4 Mq yr 1 (red solid line). The fluxes are continuum normalized. 
The main spectral lines are identified, with green and purple ticks marking the laboratory wavelength of H and Fe ii transitions, respectively. Details on the 
modelling approach can be found in G12a. We note that a stellar core radius of 300 Rq was used for the 2.5 x 1 0 4 Mq yr model (versus 60 Rq when 
M rlA = 8.5 x 1 0 4 Mq yr~* ) in order to keep the He i 7.6678 line contribution from i; A within the observed levels. 


thickness and density of the shell will cause it to break apart quicker 
after periastron as well. The observed duration during which any 
increased absorption is detected in LOS around periastron should 
thus also decrease with M IIA . 

4.4.2 Comparison to recent observed changes 

Mehneret al. (2011, 2012) show that between ~ 1998 and 2011, the 
strengths of various broad stellar- wind emission lines (e.g. Ha, H<5 
and blends of Fen, [Feu], Cru, [Cr n] ) for our direct view of the 
central source of i] Car decreased by factors of 1.5-3 relative to the 
continuum, while at the same time the line strengths for reflected 
spectra corresponding to an ;; A stellar latitude of ~75° showed no 
major variations (figs 6 and 7 of Mehner et al. 2012). The terminal 
velocity was found to be similar at all stellar latitudes and relatively 
unchanged as well (Mehner et al. 2011, 2012). The lack of any 
significant changes in )y Car’s reflected polar spectra poses a serious 
problem for the hypothesis that there has been a recent large global 
modification of M nA . Spherically symmetric ID radiative transfer 
models of ;; Car should be representative of the stellar spectrum 
at high latitudes, which are mostly unaffected by the WWC cavity 


and radiation field (G12b). To illustrate expected changes, we 
computed a ID cmfgen model similar to those described in H01, 
H06, G12a, but with decreased to 2.5 x 10~ 4 MQyr~' and 
R* iI)a increased to 300 Rq. This model shows that a factor of 3^1 
decrease in M, 1A should produce obvious changes in the reflected 
polar spectra (Fig. 12). The H<5 emission should decrease by a factor 
of ~2.5, while Ha and H/j should show a decrease in emission of a 
factor of 3-4. The broad Fe u/Crn blend between 4570 and 4600 A 
and other Fen lines should practically vanish. The emission and 
absorption of numerous other lines, like He i, would also be affected. 
Based on these results and the observed lack of any changes over the 
stellar poles, any recent decrease in (i/, A was probably by a factor 
of <2. It is beyond the scope of this paper to provide a detailed fit 
to the post-2009 observations of ij A , since other stellar parameters 
could have changed in addition to M, 1A . Such detailed analysis 
will be reported elsewhere (Groh et al., in preparation). Additional 
evidence for the lack of a recent large drop in A/,, a comes from 
/7ST/STIS measurements of the Ha EW at Weigelt blobs C and D, 
which show a long-term decrease in emission strength of only ~ 10- 
20 percent (Gull et al. 2009; Mehner et al. 2012). The continuum 
flux at Weigelt D also has not greatly changed, even though in LOS 
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Figure 13. Slices in the orbital plane from an M rjA = 2.5 x 1 O' 4 Mq yr 1 3D SPH simulation of r] Car showing the density structure (left-hand panel) and 
photoionization region created by ijb (right two panels) at <j> = 3.07, assuming H° ionizing fluxes of 10 49 ' 48 photons s -1 (middle panel) and 10 48 ' 2 photons s - 
(right-hand panel). The method for computing the photoionization region is the same as that in M12, to which we refer the reader for details. 


to r] A it has increased by a factor of 10 since 1998 (Gull et al. 2009; 
Mehner et al. 2012). 

4.5 Effects in the UV 

G12a showed that the WWC cavity in r) Car also has a big influence 
on the observed UV spectrum, which is strongly dominated by 
bound-bound transitions of Fe n. Spherically symmetric radiative 
transfer models have difficulties fitting the observed UV spectrum, 
even at apastron (H06; G12a). However, the presence of the WWC 
cavity, introduced in the context of a 2D model, reproduces well 
both the observed UV and optical spectra of i; Car by reducing the 
optical depth of Fe u transitions in LOS to i; A (G12a). 

Given the extremely large UV-emitting region of i) A (H06), any 
large change in the half-opening angle of the WWC cavity (of >20°, 
G12a) will affect the observed UV spectrum. This is in addition to 
any other changes in the UV caused by a decrease in M, 1A . Fig. 1 1 of 
H06 shows that a factor of 4 drop in A/, a drastically reduces both the 
emission and absorption of lines between 1250 and 1450 A. Note 
also the dramatic increase in emission strength (nearly a factor of 
4) of the lines between 1390 and 1410 A. Detailed 2D models of 
)) Car’s UV spectrum for different M, ]A and taking into account 
the WWC cavity are needed to better determine what changes are 
expected and can be observed. 

Mehner et al. (2011) note that the depth of the minimum in the 
UV during the 2009 photometric event, as observed in the F250W 
and F3501V filters on HST, was approximately twice as deep as 
the minimum of the 2003.5 event (their fig. 1), speculating that the 
increased depth of the 2009 event may be due to a sudden increase 
in ); A ’s mass-loss during periastron. However, it is unclear what 
would cause such enhanced mass-loss. One would normally expect 
a decrease in M J]A to reduce the size of the UV photosphere and 
make it hotter (Mehner et al. 2011). The hypothesized increase in 
mass-loss around periastron may instead be due to an increase in 
density in LOS caused by the wind of ?/ A that flows towards the 
observer during periastron passage. The primary wind between ?j A 
and the observer at periastron may now be blocking a larger relative 


amount of flux in LOS, leading to an increase in the depth of the UV 
minimum. The increase in WWC opening angle will, however, also 
alter greatly the amount of UV flux that can escape from the hotter, 
inner regions of r\ A s wind (Madura & Groh 2012). It is unclear how 
the UV flux in LOS would be affected by such a change, especially 
given the unusual circumstellar extinction. Therefore, it is difficult to 
determine whether a decrease in M, lk is consistent with the observed 
increased minimum-depth in the UV photometry. Nevertheless, the 
results of Madura & Groh (2012) show that any large change in 
WWC opening angle should affect the time-dependent directional 
UV illumination of circumbinary material in the system (see also 
Smith et al. 2004b). 

4.6 The spatially extended forbidden line emission 

Gull et al. (2009) presented an analysis of HST/STIS observations 
taken between 1998 and 2004, identifying spatially extended (up to 
0.8 arcsec), velocity-resolved forbidden emission lines from low- 
and high-ionization J species. Using 3D SPH simulations, M12 
developed a 3D dynamical model for the broad [Fein] emission 
described in Gull et al. (2009) and showed that the broad high- 
ionization forbidden emission arises in portions of the dense, ex- 
tended WWC structures photoionized by ;; B - The first detailed spa- 
tial maps showing the time evolution of the broad [Fe hi] emission 
(Gull et al. 201 1) provide strong support for the results of M12. 

M12 found that the observed broad forbidden line emission de- 
pends strongly on M, IA and the ionizing flux from r/ B . If the flux 
from ?)b remains constant, but M rA drops by a factor of 2 or more, 
the photoionization region created by should grow considerably 
(Fig. 13). The spatial extent and flux of the observed high-ionization 
forbidden lines should thus change. The emissivity should decrease 
with lii, 1A since, for densities much less than the line’s critical den- 
sity, the emissivity scales with the density squared (M12). The phase 

3 Low and high ionization refer here to atomic species with ionization po- 
tentials (IPs) below and above the IP of hydrogen. 13.6 eV. 
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dependence should differ as well, with the high-ionization emission 
vanishing at later phases (compared to earlier orbital cycles) when 
going into periastron, and reappearing at earlier phases afterward. 
Changes in the WWC opening angle will result in the broad for- 
bidden line emitting regions exhibiting different Doppler velocities, 
leading to observed variations with slit position angle and spatial 
location on the sky. New components at different spatial locations 
and Doppler velocities may appear, while components identified in 
earlier observations might vanish. 

The ring-like [Fe m] emission structures identified by Gull et al. 
(2009) and modelled by Ml 2 should disappear in future observa- 
tions if M I)A drops by a factor of 2 or more. This is because the 
dense arcs of r/ A wind where the broad high-ionization forbidden 
lines form disperse sooner after periastron the lower the value of 
M, )a (see Section 3.2). The flux, location and velocity of broad, red- 
shifted, lower-ionization lines of [Fe ii] that form in the periastron 
side of r) a’s wind (Gull et al. 2009; Teodoro et al. 2013) should 
also change if M rlA drops. The flux should decrease due to the lower 
wind density, while the spatial distance from r) A and overall out- 
ward velocity of the emitting material should increase due to the 
increased momentum of 773 ’s wind relative to that of r] A . 

During periastron passage, the high-ionization forbidden line 
emission vanishes, both the broad features in the interacting stellar 
winds and the narrow features in the Weigelt blobs (Vemer et al. 
2005; Gull et al. 2009; Mehner et al. 2010). Gull et al. (2009) found 
that the high-ionization forbidden emission was absent during the 
2003.5 event between ip = 1.0 and 1.068, but reappeared strongly 
by <p> = 1.1 22 . 4 Thus, there is an effective shutting-off of the UV 
flux from fya as it becomes temporarily enshrouded in the dense, ex- 
tended wind of j 7 A around periastron. Therefore, during the 2003.5 
event, M nA was large enough to significantly reduce the size of the 
photoionization region created by ? 7 b , from a few thousand au across 
to just several au wide. An M r]A ~ 10 ~ 3 Mg yr~' is consistent with 
this result, assuming rj B is an 05 giant with T t g *=» 40 000 K and 
an H° ionizing flux of 10 494S photons s _1 (Martins, Schaerer & 
Hillier 2005; Verner et al. 2005; Madura 2010; Mehner et al. 2010, 
M12). A factor of 4 drop in M, 1A would create a much larger pho- 
toionization region during periastron passage (e.g. middle panel of 
Fig. 13), resulting in large amounts of high-ionization forbidden 
line emission. 

In order to avoid a large photoionization region at <p w 1.068 
with such a low M, A , one must drastically reduce the number of 
ionizing photons from by a factor of ~20 (right-hand panel of 
Fig. 13). This has important implications for the spectral type of ? 7 b. 
According to Martins et al. (2005), there are no O-type giants or 
supergiants with hydrogen ionizing fluxes near 10 432 photons s _1 . 
Rather, an 08.5-09 main-sequence star would be required, with a 
T eff ss 32 000-34 000 K and log L/Lq « 4.75-4.85 (Martins et al. 
2005). However, such a star lies outside the range of acceptable 
candidates for 77 b from both the works of Verner et al. (2005) and 
Mehner et al. (2010) (see their fig. 12). Thus, a large drop in M, 1A 
should reveal itself in the broad forbidden line emission during 77 
Car’s next periastron passage in 2014. If not, the implication is 
either (1) M^ A has not dropped significantly since 2004 or (2) our 
current understanding of r/B is incorrect and the star has a lower 
effective temperature and hydrogen ionizing flux. 

Another major issue regarding the forbidden line emission con- 
cerns the population of available Fe + to be photoionized to Fe 2+ 

4 These are spectroscopic phases, which are not necessarily equal to the 
orbital phase. The two are not expected to differ by more than four-to-six 
weeks. 


by 77 B so that [Fein] emission can form. It is generally accepted 
that the energy required to ionize Fe + to Fe z+ comes from photons 
provided by rj B (Verner et al. 2005; Gull et al. 2009; Mehner et al. 
2010). Thus, the iron in any [Fe m] emitting region needs to initially 
be in the Fe + state. H01, H06 show that the ionization structure 
of tja’s wind and the amount of Fe + present depend strongly on 
charge exchange processes between Fe 2+ and neutral H via the re- 
action Fe 2+ + H -o- Fe + + H + . According to fig. 9 of H01, which 
assumes M nA = 10~ 3 Mg yr _1 , Fe + is the dominant ionization state 
of iron for radii r > 50 au in tj A ’ s wind. Below this, higher ioniza- 
tion states dominate. However, in cmfgen models of r/ A that assume 
M nA = 2.5 x 10~ 4 Mq yr~' (with the other parameters held fixed), 
H recombination never occurs and the iron is primarily Fe 2+ in the 
outer wind. H01, H06 discuss this issue and use it as a reason for 
favouring a higher M, IA , as it provides a better fit to the observed 
spectrum. 

Thus, if M, 1a were to drop to ~2.5 x 10 ~ 4 Mq yr~' , Fe 2+ would 
be the dominant ionization state of iron in the outer 17 A wind. f/ B 
would then be unnecessary for [Fe 111 ] emission to form, provided 
the 77 a wind is in the right density and temperature regimes and has 
not recombined to Fe+ or been excited to higher states. This has a 
number of major ramifications. First, while a hot 773 would not be 
necessary for [Fe 111 ] emission, 77 s ’s UV flux would be necessary for 
forming the other observed higher ionization lines of [Ne m] and 
[Arm] (Vemer et al. 2005; Gull et al. 2009). A large drop in M r/A 
would hence lead to huge differences between the wind structures 
observed in lines like [Arm] and those of [Fein], 

Secondly, if Fe z+ became the dominant state of iron, then broad, 
red-shifted [Fe 111 ] emission would become observable in most di- 
rections for most of the binary orbit, arising in the periastron side 
of 77 A ’s wind and portions of the WWC region that are moving 
away from the observer. More importantly, the broad [Fe m] emis- 
sion would not completely vanish during periastron passage. The 
sudden appearance of large amounts of [Fe 11 ] emission during pe- 
riastron passage, as previously observed (Gull et al. 2009), would 
also not occur. If the photoionization region extended far enough, 
the narrow high-ionization lines that arise in the Weigelt blobs and 
which also vanish during periastron would remain visible as well. 

It is conceivable that in such a low M, 1A situation, the Fe 2+ in 
? 7 a wind that gets compressed in the WWC regions will recombine 
to Fe + . Such material on the 77 a side of the system could then be 
photoionized to reform Fe 2+ and allow [Fe m] emission. Iron in the 
dense shells of compressed primary wind on the periastron side of 
?7a created during periastron passage could also recombine to Fe + . 
As these shells would be shielded from 773 ’s radiation, the gas in 
them would remain in the Fe + state and produce broad red-shifted 
[Fe 11 ] emission. However, the uncompressed 77 a wind would still 
consist of Fe 2+ and exhibit significant amounts of broad red-shifted 
[Fe 111 ] emission. This emission would appear at nearly all phases 
in directions to the south-east on the sky at Doppler velocities up 
to Uoo.tja ■ The latest HST observations show no such emission (Gull 
et al. 2011). However, future HST observations of the forbidden 
lines could help constrain any new changes to M, 1A . A major benefit 
of studying the extended, broad forbidden line structures is that they 
provide an excellent record of 77 Car’s mass-loss over the past 2-3 
orbital cycles (Gull et al. 2009, 2011; M12; Teodoro et al. 2013). 

4.7 Other expected changes 

There are other observable changes expected from a large decrease 
in M JtA that are worth briefly mentioning. First, a drop in M 1]A could 
affect any periodic dust formation during periastron passage, as 
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mentioned in Smith (2010). A large decrease in M nA would decrease 
the density in the compressed WWC zone where dust is expected 
to form, leading to a decrease in the near-IR flux and/or changes 
to the spectral energy distribution. An overall decrease in primary 
wind density and WWC-region density could also allow the intense 
radiation from the stars to more easily destroy any dust that forms. 
Older dust that lies in extended regions farther from the stars could 
also be destroyed or heated if intense radiation flowing through the 
WWC cavity can reach it. This may help explain the recent decrease 
in extinction and increase in brightness of the central source in LOS, 
e.g. an increase in WWC opening angle that allows more radiation 
to escape in LOS and destroy dust between us and the stars, thereby 
altering the extinction in LOS. 

The wavelength-dependent apparent size of r) A is also controlled 
by the wavelength- and radially dependent opacity and emissivity, 
which depend on M, A (HOI, H06, G12a). As M nA drops, the contin- 
uum extension of 17 A in the optical and near-IR should drop as well. 
In other words, tj A should appear smaller than it was when observed 
between 2003 and 2004 (van Boekel et al. 2003; Weigelt et al. 2007; 
Groh et al. 2010b). Continuum flux changes at a given wavelength 
that occur when the apex of the WWC cavity penetrates the surface 
of optical depth r = 2/3 (e.g. the ‘bore-hole’ effect, see Madura 
2010; Madura & Owocki 2010) would also be affected, changing 
the behaviour of the optical and near-IR photometric eclipse-like 
events observed to occur every 5.54 yr around periastron (White- 
lock et al. 2004; Femandez-Lajus et al. 2010; Madura 2010). A 
reduced M VA should cause the eclipse-like events to start later due 
to the longer time it takes for the WWC cavity to reach optical 
depth r = 2/3 in the wind photosphere of ;; A - Similarly, the events 
should end earlier as the apex recedes from the location of optical 
depth r = 2/3. For low enough M,, A (< 10 -4 Mq yr -1 ), the WWC 
cavity would never penetrate ;j A ’s optically thick wind photosphere 
and the photometric eclipse-like events would very likely cease. As 
of the writing of this paper, no significant decrease in the apparent 
continuum extension of r) A in the near-IR has been observed. No 
major differences are seen between the 2003.5 and 2009 events in 
the ground-based optical photometry either (Femandez-Lajus et al. 
2010 ). 


5 SUMMARY AND CONCLUSIONS 

Using 3D SPH simulations, we showed that M,, A has profound 
effects on the 3D time-dependent hydrodynamics of 17 Car’s binary 
colliding winds on a wide range of spatial scales. These simulations 
lay the foundation for future work generating synthetic data that 
can be compared directly to past and future observations of r] Car. 
Such modelling provides information that is important for helping 
constrain r] Car’s recent mass-loss history and possible future state. 

The 3D SPH simulations and ID cmfgen models of ?; A ’s spectra 
show that for the range of explored values, a factor of 2 or more 
drop in M r]A results in substantial changes to numerous observables 
across a wide range of wavelengths. Below we summarize the most 
important changes expected. 

(1) For A/, a ~ 10 -3 Mq yr -1 , our LOS intersects portions of the 
dense, outer WWC region. For lower M, A (by factors of 2 or more), 
our LOS passes primarily through the wind cavity carved by 
An increase in WWC opening angle may thus explain the observed 
brightening of the system and changes to various wind-emission 
line profiles seen in LOS, with little to no changes at high stellar 
latitudes or at the Weigelt blobs. 


(2) Increases in WWC opening angle cause LOSs to the WWC 
apex to move into r] A s dense wind at later 0 as the system goes into 
periastron, and to move out of ;; A ’s wind at earlier cj> as the system 
leaves periastron. Combined with the decreased cooling-transition 
phase of r/ b’s post-shock wind, the net result is a shorter X-ray 
minimum (by ~40-80 d) during periastron passage for simulation 
Cases B and C, versus Case A. 

(3) A decrease in M 11A leads to a corresponding drop in X-ray flux 
at phases far from periastron. Decreases in M rlA also propagate into 
fj Car’s extended interacting winds and may change the thermaliza- 
tion and/or absorption of the CCE component. 

(4) Increases in WWC opening angle cause the trailing arm of 
hot post-shock 77 b gas to cross the observer’s LOS at later (f> as the 
system approaches periastron, delaying when the strongest observed 
X-ray lines become increasingly blueshifted and broader. 

(5) As M VA is lowered, the spatial extent of the He + zone in 
r] A s dense inner wind increases considerably, altering the timing 
of the emission peaks observed in He ii 7.4686 around periastron. 
Assuming an initial M VA = 10 -3 Mq yr -1 , a factor of 4 drop should 
lead to the detection of significant He n emission at all orbital phases. 
A factor of 2 drop should change the timing of the He ii event, with 
the first peak occurring ~3 months earlier and the second peak 
lasting ~2-3 months longer. 

(6) An increase in WWC opening angle will greatly affect 
broad wind-emission features of Ha, H/3 and Fen lines. The 
amount of Balmer P-Cygni absorption in LOS will be reduced 
and possibly entirely absent. The latitudinal and azimuthal vari- 
ation of the Balmer P-Cygni absorption will also change. For 
M r]A < 2.5 x 10 -4 Mq yr -1 , there could be little to no absorption 
in reflected spectra over the stellar poles of tj A . 

(7) The observed increase in Balmer P-Cygni absorption in LOS 
during 17 Car’s spectroscopic events should decrease for decreased 
M, a . The density and amount of primary wind that flows to- 
wards the observer during periastron passage will also decrease. 
For M Va < 2.5 x 10 -4 Mq yr -1 , the absorption in LOS during the 
events should vanish since 77 A will be able to easily ionize hydrogen 
throughout its extended wind. 

(8) A decrease in M r]A by a factor of 2-4 should decrease the 
emission strength of Ha, H/3 and H<5 by a factor of 2-3. The broad 
Fe n/Cr ii blend between 4570 and 4600 A and other Fe ii lines should 
practically vanish. Such changes should be visible not only in LOS, 
but also at high stellar latitudes and in other spatial directions, such 
as at the Weigelt blobs. 

(9) Given the extremely large UV-emitting region of r/ A , large 
changes in WWC half-opening angle of ~20° will alter the observed 
UV spectrum. Detailed 2D radiative transfer models of the resulting 
spectrum are needed, but current ID models show that a factor of 
4 drop in M J]A reduces both the emission and absorption of lines 
between 1250 and 1450 A, and increases the emission strength of 
lines between 1390 and 1410 A. 

(10) The observed broad forbidden line emission depends 
strongly on M nA and the ionizing flux from tie- As M iia drops, 
the size of the photoionization region where the highest ionization 
forbidden lines form grows considerably. The size, location, flux 
and time variability of the various forbidden lines will thus change. 
Ring-like emission structures that form in the dense arcs of pri- 
mary wind on the apastron side of the system will be less structured 
and could vanish altogether since the arcs disintegrate sooner after 
periastron for lower Af, A . 

(IT) A decreased 77 A wind density due to a large drop in M J]A may 
lead to the inability of ij A s wind to trap the ionizing radiation from 
77b at periastron. This would cause a drastic change in the behaviour 
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of the low- and high-ionization forbidden emission lines during pe- 
riastron passage. The high-ionization lines would no longer vanish 
at periastron, while the low-ionization lines would fail to appear. 
This behaviour could occur in both the extended wind structures of 
tj Car and the Weigelt blobs. 

(12) For M 1a < 2.5 x lCL 4 MQyr~ 1 , Fe 2+ becomes the domi- 
nant ionization state of iron at large radii in ?? A ’s wind, resulting in 
large amounts of broad, red-shifted [Fe in] emission in most direc- 
tions over most of the orbital cycle. 

(13) A large decrease in M 1JA could affect periodic dust formation 
during periastron passage since it would decrease the density in 
the compressed WWC zone and increase the likelihood of heating 
and/or destruction by the intense stellar radiation fields. 

(14) The wavelength-dependent apparent size of a in the op- 
tical and near-IR continuum decreases with M, )A . Continuum flux 
changes that occur when the WWC apex penetrates the surface of 
high optical depth in ;y A ’s wind photosphere would be affected, 
making the periodic photometric eclipse-like events shorter due to 
the smaller continuum size of ij A . For M nA < 10~ 4 Mq yr _1 , the 
eclipse-like events would likely cease. 

Given all of the multiwavelength changes expected, a recent fac- 
tor of 2 or more decrease in M rlA should be discernable with the 
right observations. However, as discussed in the previous section, 
there are a large number of inconsistencies between what is ob- 
served and what numerous models predict should occur if M, 1A 
were to drop by a factor of 2^1 from its historically accepted value 
of ~1CL 3 Mq yr~' . One of the most problematic is the lack of any 
observed changes in the broad wind-emission lines at high stellar 
latitudes or different spatial locations (Mehner et al. 2012). Also 
puzzling is the fact that the timing of the start of the 2009 spectro- 
scopic event in X-rays, He n L4686 and the ground-based optical 
photometry was not delayed by any noticeable amount. It further- 
more seems odd that M nA would change with no corresponding 
change in wind terminal speed. Taken together with the other un- 
observed but expected changes, it appears that contrary to earlier 
claims, there has not been a recent significant drop in A/, a . Current 
evidence suggests that any decrease in M J]A was likely by a factor 
of <2 and occurred sometime after 2004, quite possibly during the 
2009 event itself. 

We speculate that most of the observed recent changes in tj Car 
are due to a slight (<10°) increase in the half-opening angle of 
the WWC zone, which altered the amount of extinction in LOS to 
the central source. This would explain why the continuum flux and 
broad wind-emission features in LOS have changed, while features 
at high latitudes and in other directions (e.g. at the Weigelt blobs) 
remain relatively constant. It would also explain why the start of 
the 2009 event at multiple wavelengths did not change. A modest 
decrease in M nA by a factor of <2 sometime after 2004 could be 
responsible, and would hardly be surprising given ij Car’s nature. 
However, changes in ii B ’s wind/stellar parameters, while less likely, 
cannot be fully ruled out at this time. Other less likely alternatives 
to a drop in are discussed in Mehner et al. (2010, 2012). 

Finally, we note that the results of Sections 3.1.1 and 3.1.2 depend 
on rather complicated shock and radiative-wind-driving physics. 
In addition to the known effects of RI and braking, there is the 
rather new phenomenon of so-called self-regulated shocks (SRSs), 
in which ionizing X-rays from the WWC shocks inhibit the wind ac- 
celeration of one or both stars, leading to lower pre-shock velocities 
and lower shocked plasma temperatures (Parkin & Sim 2013). For 
most of r] Car’s orbit, the SRS effect is likely to be negligible due 
to the large binary separation. However, during periastron passage, 


SRSs may further inhibit the acceleration of tj B ’s wind, helping 
trigger a cooling-transition phase of tjs’s post-shock wind and a 
possible WWC ‘collapse’. SRSs may cause the cooling-transition 
phase to occur earlier and/or last longer compared to what our 
simulations predict. The potential for a WWC collapse on to ijb 
in the context of SRSs is difficult to determine without a detailed 
time-dependent model since radiative braking (not included in the 
analysis of Parkin & Sim 2013) may prevail over the SRS effect 
and prevent a wind-photosphere collision (Parkin & Sim 2013). 
Inclusion of SRSs in future models of r) Car is thus worthwhile and 
may further elucidate the peculiar behaviour of X-rays, He ii 7.4686 
and other lines around periastron. 

Many recent works suggest that very massive stars experience ex- 
treme r) Car-like mass ejections of ~ 1-10 Mq in the year to decade 
before they explode as Type II SNe (Kotak & Vink 2006; Smith 
et al. 2007; Ofek et al. 2007, 2010, 2013; Chevalier 2012; Quataert 
& Shiode 2012). The narrow emission lines indicative of Type Iln 
SNe are interpreted as signatures of the SN ejecta interacting with 
this freshly expelled circumstellar material (Smith et al. 2007; Ofek 
et al. 2007, 2013). In fact, the most likely detection of a progeni- 
tor before its explosion as a Type Iln involved an LBV (Gal- Yam 
& Leonard 2009; Ofek et al. 2013), thus raising the question of 
ri Car’s near-term fate, Car’s numerous massive eruptions point 
strongly to repeated phases of instability and may be an indication 
of its imminent demise as a powerful Type Iln SN (Smith et al. 
2007; Stritzinger et al. 2012). Constraining Car’s recent mass- 
loss history is therefore crucial for determining its possible near- 
and long-term states. The 3D simulations and analysis in this paper 
provide detailed information that is important for helping achieve 
this goal. Extensive monitoring of 17 Car’s next periastron passage 
in 2014 will undoubtedly play a pivotal role as well. 
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APPENDIX A: AN ANTIGRAVITY APPROACH 
FOR IMPLEMENTING RADIATIVE FORCES 
IN SPH SIMULATIONS OF COLLIDING 
STELLAR WINDS 

Al Basic approach 

The SPH method used for the simulations in this paper follows the 
Lagrangian motion of the modelled gas particles, removing the need 
for a complex spatial grid. In the SPH scheme, spatial derivatives 
are calculated via analytical differentiation of interpolation formu- 
lae (Monaghan 1992), as opposed to Eulerian grid-based methods 
where derivatives are computed using finite-difference or finite- 
volume techniques. While the SPH approach allows us to perform 
large-scale 3D simulations with significantly less computational 
cost compared to a grid-based code, the lack of a spatial grid com- 
plicates the calculation of the projected velocity gradient required 
for computing the CAK (Castor, Abbott & Klein 1975) line force 
that drives the individual stellar winds. We therefore employ an 
alternative approach, described below, that produces a stellar wind 
that follows the desired beta- velocity law and is consistent with the 
results of standard CAK theory. 
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We start by considering the radial equation of motion for a steady, 
spherically symmetric stellar wind driven by a radiation force g ra d = 
k L t /(4ncr 2 ) and competing against stellar gravity GM+fr 2 '. 


Av 

V— = grad - 
dr 


GM* 


(Al) 


Here, v is the wind speed, r is the local radius, M t and L* are 
the stellar mass and luminosity, k is the wind opacity and c is 
the speed of light. We neglect standard sound speed terms that 
are, in practice, negligible. Since the key to driving a stellar wind 
is overcoming gravity, it is convenient to define a dimensionless 
equation of motion that is scaled by the gravity, 


du> 

dx 


1 , 


(A2) 


where the dependent variable w = v 2 /vl sc is the ratio of the ra- 
dial kinetic energy to the surface escape energy v 2 sc = 2 GM t / R„, 
and the independent variable is the inverse radius coordinate x 
= 1 — R t /r, measured from the stellar surface radius R„. T = 
dgiadKGM*) is the ratio of the local radiative force to gravity, 
sometimes called the ‘Eddington parameter’. 

A simple model to consider is one in which we assume that, at 
some R„, g rad suddenly jumps from below gravity to some fixed fac- 
tor T > 1 above gravity. In such a simple ‘antigravity’ model, static 
material at the surface R, effectively ‘falls’ away from the star with a 
local net outward acceleration at radius r of (r — ljGM^/r 2 . In terms 
of equation (A2), integration yields a solution w(x) = (r — l)x, 
which in dimensional form gives u(r) = i>oo(l — R t /r) 1 ^ 2 , where 
Utxj = UescVr — 1. This is just a specific p = 1/2 form of the beta- 
velocity law v(r) — Voo(l — R t j r'f . 

While we considered only a simple example above, this anti- 
gravity approach can readily give a beta- velocity law with any index 
p. To demonstrate this, we rewrite equation (A2) in terms of v(x) 
and require the Eddington parameter to have a spatial variation T (x). 
Then, 


2v(x) d 
dx 


v(x) = T(x) - 1. 


(A3) 


With v{x) = VqqX^ and defining = u^/tr, c , the above gives 
T(x) = 2 (Suwr^ 1 + 1 . (A4) 


we require that the ratio of the two stellar luminosities satisfy 



(A8) 


where a, typically = 2/3, is the CAK power index. Given an ob- 
served primary star luminosity L \ , the above equation can be used 
to find that of the secondary, L 2 . Therefore, with L\ and M,, if,, Mi, 
/I, and , specified for each of the stars, the total force per mass 
g tot exerted on the wind particles can be computed as corrections to 
the vector gravities g, of the stars, 


g lol 



HLj \ 
\nGM\c ) 


gi + 



icL 2 \ 
4nGM 2 c ) 


gi- 


(A9) 


With the above formalism in place, we incorporate the effect of 
radiative forces into the 3D SPH simulations by 


(1) Specifying the input parameters L\ = 5 x 10 6 Lq and M,, 
if,, Mj, /3i and , for each star (Table l). 

(2) For the given L\ and assumed stellar mass and mass-loss rate 
ratios, using equation (A8) with a = 2/3 to obtain L 2 . 

(3) Using equation (A5) to compute K\ and k 2 . Particles from 
each star are distinguished via their specified distinct masses. 

(4) Computing k at each location in the wind using equation (A6). 

(5) Computing the associated g to , at each point with equa- 
tion (A9). 


The stellar winds are initiated at twice the sound speed in order 
to ensure information characteristics that point from the star into 
the computational domain. The SPH particles from each star i are 
injected at a radius finj,,- = if, / [ 1 — (2c s /v 0Ot i) l ^ i ], where the sound 
speed c s = -JT yk^/ p and T, y, k% and /r are the photospheric 
temperature, ratio of specific heats, Boltzmann’s constant and mean 
molecular weight, respectively. 

Fig. Al shows that the antigravity method produces the cor- 
rect wind-velocity and density structures. Plots of v(r) and p(r) 
through the back side of ;?a’s wind, taken at apastron from the 
M, a = 8.5 x 1(T 4 Mq yr~' 3D SPH simulation, are compared 
to the analytic expressions v(r) = Doo(l — R+/r)P and p(r) = 
M /\\nr 2 v{r)\ for the wind structure of an isolated star, assum- 
ing the values for r ] a from Table 1. The SPH simulation results 
are in very good agreement with the analytic scalings, which helps 


Obtaining the desired velocity laws simply requires specifying the 
correct opacity for each wind. Using equation (A4) together with 
T(r) = r 2 g rad /(GM*) = /cL*/(47tGM*c) yields 

4-7TC 

Ktirt) = — [<,Tf/A(l - Ri/r tf*- 1 + GMi] , (A5) 

where i — 1,2 correspond to the individual stars and r, is measured 
with respect to star i. Achieving the different wind terminal speeds 
requires K\ ^ k 2 . However, regions where the two winds mix re- 
quires weighting the local opacity by the relative mass contribution 
from each star. 


k =Ki +k 2 = K i +(K2-Ki)^ . (A6) 

where p = p\ + p 2 is the total density. 

In order to ensure that the stellar mass-loss rates correspond to 
what is expected from CAK theory, e.g. 


4/cak = 


L t a 
c 2 1 — a 


Q r 


- (1— a)/a 


l - r 


(A7) 
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Figure Al. Line plots of the total wind speed (orange and black) and 
log density (red and blue) as a function of distance from rj,\ (in units of 
a = 15.45 au) taken from the M, A = 8.5 x 10~ 4 Mq yr^ 1 3D SPH simu- 
lation at apastron (solid lines) compared to the simple analytic expressions 
for v(r) and p(r) for an isolated star (dashed lines). 


(Gayleyl995), 
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Figure A2. Line plots of the opacity <e; (equation A5) as a function of 
distance from the star (in units of a = 15.45 au) for rj A (top) and t; a (bottom) 
for the M, )a = 8.5 x 10 4 Mq yr _l simulation. 


verify our approach and results. Fig. A2 shows line plots of the 
opacity at, used to drive the individual stellar winds for the desired 
mass-loss rates and wind terminal speeds. Note that k 2 is signif- 
icantly larger than K \ . More extensive testing and analysis of the 
antigravity method can be found in Russell (2013). 


A2 An analytic solution for the flow along the line of centres 

Since the radiative force from both stars is present throughout the 
computational domain, the velocity-altering effects of RI and stellar 
gravity are included in the simulations (Stevens & Pollock 1994; 
Russell 2013). To assist in the interpretation and understanding 
of the 3D SPH simulations, we derive here a sample analytic ID 
solution for the effect of the antigravity radiative forces on the 
flow along the line of centres between the stars. Because the SPH 
simulations use /S = 1 velocity laws to parametrize the stellar winds, 
we focus on the /) = 1 case. 

As yib has a much lower luminosity than q A , its wind will be the 
one most strongly affected by the radiative force of its companion. 
We therefore investigate how q b’s wind velocity along the line of 
centres is affected by the radiation and gravity of q A . The equation 
of motion along the line of centres in terms of the radius r measured 


from jjb , for a given binary separation D , is 


du 2 GM\ 

Vi = - 

dr (D — r ) 2 


1 - 


kL\ 


AnGM\c 


G Mi 


1 - 


tcL 2 \ 
AnGMic ) 
(A 10) 


where subscripts 1 and 2 refer to q A and ?7 b , respectively. Because 
we are interested in the effect of q A in slowing (/b’s wind before the 
WWC, we can set ic — k 2 - Multiplying both sides of equation (A10) 
by r 2 /(GMi) and using the definition of the Eddington parameter 
F; = KjLj /(4nGMjc), we get 


v 2 r~ dv 2 

GM 2 dr 


= (F 2 — 1) — 


(D - r) 2 


Mi 

~Mi 


k 2 T i 
K 1 


- 1 


(All) 


Defining x = 1 — R 2 /r, w = {v\R 2 )/(/2GM 2 ), q = L\/L 2 , m = 
Mi/M 2 , and using the fact that (c,- = (4rcG A/; cT, ■)//,;, equation 
(A1 1) becomes 


dui 

— = (r 2 -i)- 

dr 


qT 2 — m 


[D(l — x)/R 2 — l] 2 


(A 12) 


Equation (A4) gives the spatial variation of T with x. Therefore, for 
/S = 1 and u>oo = (tv 2 /tV, 2 ) 2 . r 2 (x) = 1 + Iw^x and 


du) g(l + 21 UOO.V) — m 

— = 2w<xX . 

dr [D(l — x)/R 2 — l] 2 


(A13) 


Equation (A13) can be integrated straightforwardly to give the 
variation of the scaled wind kinetic energy w as a function of the 
scaled inverse radius x. Letting d = D/R 2 , we find 

IqWoo ln[(l — x)d — 1J 


U){x) = WqoX — 


d 2 


qd — md — 2qWoo + 2 qdw a 
d 2 [\ +(x - 1 )d] 


(A 14) 


With w/Woo = (Vi/VaO'i) 2 , the radial variation of the velocity of 
? 7 b "s wind along the line of centres can be found from 


= 1 


u oo,2 


Ri 

r 


+ 


2qR\ 
" D 2 

2 GM 2 (q — m) 

vlo,2 D C ~ D /C 


In 


D 


- 1 


2qR 2 (\ - f ) 
D(1 - D/r) 


(A15) 


The corresponding expression for the radial variation of q A s wind 
velocity iq along the line of centres is similarly 


■Vi 


fl-^V-^Mn 

1 

1 

Q 1 

\ r' J D- 

V r ’ 


+ 


2q'R l (l - f ) 
D(1 - D/r') 


2 GM\{q' — in’) 

+ vl A D(l-D/r')' 


(A16) 


where q' = 1 /q, in' = 1 /in and r 1 is measured from q A . 

In addition to the effects of RI, orbital motion will increase the 
wind speeds as the system moves from apastron to periastron. The 
stellar velocity along the line of centres for each star i with respect 

to the system COM is v r j = \J vl — u 2 ; , where 
v t | = h Mi/ (Mj D) 


v t , 2 = hMi/(M T D) (A17) 

are the transverse components of the stellar velocities with respect 
to the COM, and M x = M t + M 2 , h 2 = GaMj( 1 — e 2 ), G is the 
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Distance from t/„ (au) Distance from t,„ (an) 


Figure A3. Total <73 wind speed 1 ) 2 , tot = t >2 + v r 1 along the line of centres 
as a function of distance from 773 assuming the parameters in Table 1 
and = 8.5 x 10 4 Mq yr , q = 22.28, for binary separations D 
corresponding to 0 = 0.5 (apastron), 0.9, 0.96, 0.98 and 0.99. For reference, 
the solid black line is the standard f) = 1 law without RI effects. Each 
curve is truncated at the radius of ram pressure balance rp given by equation 
(A19), indicated by the vertical lines. There is no stable balance at periastron 
(0 = 1 ), implying that 773 cannot drive a wind towards tja - 


gravitational constant, a is the semimajor axis length, e is the orbital 
eccentricity and t>* is the stellar orbital speed given by 


v~ — GMj 




(A18) 


Fig. A3 plots the total 773 wind speed along the line of centres 
v 2t t ot = D 2 + tV, 2 versus r for several D corresponding to various 
phases leading up to periastron. Each curve is truncated at the radius 
of wind-wind ram pressure balance 7 - b measured from 773 (indicated 
by the vertical lines), which is found by solving the ram pressure 
balance condition 


D I Mi[vi(r' = D — r b ) + v r 1 ] 

- = 1 + a — , (A19) 

/' b y M 2 [v 2 (r = r b ) + v r , 2 ] 

where v 2 (r = r b ) and v\(r' = D — r b ) are given by equations (A15) 
and ( A 1 6) . The standard /3 = 1 velocity-law profile for 773 is included 
for comparison. 

Note how when the effects of RI are included, 773 ’s wind is 
initially accelerated from the stellar surface, but then eventually 
decelerates as it approaches The offset position from 77 A of 
the deceleration depends on the luminosity ratio q. This decel- 
eration weakens the colliding wind shocks in the 3D SPH sim- 
ulations and enhances the dominance of ij A s wind over that of 
773 , thus moving the WWC apex closer to 773 and narrowing the 
opening angle of the WWC zone around 773 . Because it has a 
strong dependence on the binary separation D, these velocity- 
altering effects are time dependent and most pronounced near 
periastron. 

We find there is no stable balance solution at 0 = 1, i.e. according 
to equation (A 15), 773 cannot drive a wind along the line of centres 
towards q A at periastron. This implies that in the absence of other 
effects, ? 7 A ’s wind would collide with the surface of 773 . However, 
because we set k = k 2 in equation (A10), there is a strong coupling 
between q A s radiation field and 773 ’s wind. To investigate the effect 
of a weaker coupling analogous to that used by PI 1, we rederived 
equation (A15) assuming k = K\ for the coupling between 7 / A ’s 
radiation and 773 ’s wind. Fig. A4 plots the resulting total 773 wind 
speed along the line of centres. One notices immediately that the 
weaker coupling decreases greatly the importance of RI and leads 
to much higher pre-shock 773 wind speeds. The higher pre-shock 


Figure A4. Same as Fig. A3, but using k = K\ for the influence of tja’s 
radiation field on 773 ’s wind and including periastron (0 = 1.0). The solid 
blue line for 0 = 0.5 overlaps almost exactly with the solid black line that 
corresponds to the standard p = 1 law without RI effects. Unlike Fig. A3, 
there is a stable ram balance at periastron. 


wind speeds reduce the importance of radiative cooling in the post- 
shock gas and there is no longer a ‘collapse’ of the WWC zone 
at periastron, but rather a stable ram balance. This is the same 
result obtained by PI 1, indicating that a strong coupling between 
77 A ’s radiation field and 773 ’s wind is necessary for RI effects to be 
important. 


A3 Possible changes to M 

Because we fix the mass-loss rates in our antigravity approach, 
possible changes to the stellar mass-loss due to RI are not included. 
Such changes though are not expected to be significant in 77 Car, 
with the exception of at periastron for ij B . 

As demonstrated in Russell (2013), the mass-loss rate factor 
7 W that describes the change in star l’s mass-loss rate due to the 
radiation and gravity of star 2 is 



1 - 


M 2 /Mi 


(D/ Ri - 1) 2 J 


(A20) 


Using the values in Table 1 and a = 2/3, we find that the radi- 
ation and gravity of 773 are insufficient to modify ij A ’s mass-loss 
rate at any phase, and so m = 1 for the entire 5.54 yr orbit and 
there is no difference in M nA compared to the expected single-star 
values. 

For most of 77 Car’s cycle (0.5 < 0 < 0.98), D is large enough 
that ?7a does not modify 773 ’s mass-loss and 777 sa 1. At 0 = 0.99, 
777 0.87 and there is a slight reduction in 773 ’s M. Only between 

0 = 0.996 (777 sa 0.66) and periastron {m sa 0.52) is there a sig- 
nificant reduction in the M of 773 . Since this occurs only for an 
extremely brief period at periastron, it should not greatly affect 
our results or conclusions. We further note that the above results 
assume M rjA = 8.5 x 10 ~ 4 Mq yr~* and q = 22.28 (Table 1). Be- 
cause our lower M, 1A simulations use slightly larger 773 luminosities, 
the reduction in 773 ’s M will be less, with m as 0.66 and 0.79 at pe- 
riastron for M, 1a — 4.8 x 10 ~ 4 Mq yr -1 and 2.4 x 10 ~ 4 Mq yr -1 , 
respectively. 
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APPENDIX B: THE STRUCTURE OF ETA 
CAR’S BINARY COLLIDING WINDS IN THE xz 
AND yz PLANES PERPENDICULAR 
TO THE ORBITAL PLANE 

BI Results from the small-domain simulations 

Figs B 1 through B4 show density and temperature slices in the xz and 
yz planes perpendicular to the orbital plane at the same five orbital 
phases of Figs 1 and 2 for the small-domain 3D SPH simulations, 
allowing one to more fully appreciate the complex 3D structure of 
the winds and their collision. Row a of Figs B 1 and B2 further 
demonstrates the axisymmetry of the shock cone at phases around 
apastron, the increase in cavity opening angle with decreasing M nA , 
and the increase in distance between the WWC apex and ?) A as M nA 
is increased. Decreases in the thickness, density and stability of the 
cold, post-shock primary wind region with decreasing M VA are also 
apparent. 

At 4> = 0.9, the WWC zone starts to become distorted and )/ B 
begins to move to the back side of r) A . Slices centred in the xz plane 
sample the leading edge of the trailing arm of the WWC zone above 
and below the orbital plane. Fig. B2 shows the changing temperature 
structure of the hot post-shock ij B wind as the trailing arm sweeps 
through the xz plane. The xz slice at <j> = 0.9 samples more of the 
leading edge of the trailing arm the higher the value of M, ]A . This 
is because higher M, IA produce smaller shock-cone opening angles. 
The overall result is the appearance that the arms of the hot shock 
increase in width in the xz plane as the system moves from <f> = 0.5 
to 0.9, filling in the cooler wind cavity. 

At periastron (row c), tj b is behind rj A , and there is a brief moment 
when the wind of r] A can flow unimpeded in the +x (apastron) 
direction. During this time, the apex of the low-density cavity on 
the +x side of the system fills with dense ri A wind. The inner r] A 
wind region in the +x direction also starts to look more spherical 
(see e.g. row d of Fig. Bl). The thickness and density of this inner 
primary wind region increases with M, 1A . The wind momentum 
balance is responsible for this since a higher value of M, 1A moves 
the apex of the WWC zone farther from ;; A , allowing the primary 
wind to fill a larger volume in the +x direction during periastron 
passage. Slices in the xz plane at periastron also sample mostly cold 
(~10 4 K) material from both winds, although the )/b wind cavity 
in the lower M nA simulations contains warmer material due to the 
larger WWC opening angle and less oblique shocks. 

Following periastron, i-] B emerges from ij A s dense wind, carving 
a new low-density wind cavity that is visible on the — x side of the 
xz-plane slices (Figs B 1 and B2, row d). Eventually r ) B moves back 
to the +x side of the system and the WWC zone is reestablished 
(row e). A new wind cavity is slowly carved in ?j A ’s wind on the +x 
side while the spiral WWC region and low-density cavity created on 
the —x side flow outward. Due to the more similar wind momenta, 
the wind cavity is larger and hotter in the lower M rA simulations. 
The size and geometry of the cavities on the +x and — x sides 
of the system are also very distinct. The wind of r] B is able to 
plough through the dense shell of primary wind on the +x side 
of the system at earlier phases after periastron in the lower M rA 
simulations (Fig. Bl, row e). 

For most of the orbit, slices of the density and temperature in 
the yz plane are very similar for all three simulations since only the 


dense wind of r] A is sampled. The exception is around periastron, 
mainly for Case C (Figs B3 and B4). There is a clear density and 
temperature difference between the +y and — y sides of the system 
in the yz-plane slice for Case C at <j> = 0.9 (row b). This is because of 
Case C’s extremely large WWC-cavity opening angle. Examination 
of row b of Fig. 1 shows that at (f> = 0.9, only in Case C has the 
leading arm of the WWC zone passed through the yz plane. In Cases 
A and B, the leading arm does not pass through the yz plane until 
close to <f) = 1 (Figs B3 and B4, row c). 

Again, following periastron passage, r) B returns to the +x side 
of the system, carving a new wind cavity along the way. We see 
in the density and temperature slices taken at <p = 1.03 (Figs B3 
and B4, row d) that this produces a difference in the size of the 
cavity on the +y and — y sides of the system. The magnitude of 
this difference increases with decreasing M, 1A . The larger cavity on 
the +y side of the system is the remnant of the WWC cavity from 
before periastron, while the smaller cavity on the — y side is part 
of the new cavity created after periastron. A temperature difference 
also exists within the cavity, with the +y side being much colder 
than the — y side in the simulations for Cases A and B. After enough 
time passes, the new cavity on the —y side of the system expands 
to a size comparable to the cavity on the +y side, although there is 
still a clear difference in the temperature and density structure of 
the two cavities (Figs B3 and B4, row e). 

B2 Results from the large-domain simulations 

FigsB5 through B8 show density and temperature slices in the xz and 
yz planes at the same five orbital phases of Figs 8 and 9 for the large- 
domain 3D SPH simulations. There is a clear left-right asymmetry 
in the density and temperature in each panel of the figures. Fig. B5 
shows that the higher the value of M nA , the smaller the wind cavities 
carved by )/ B on the — x side of the system become in the xz plane. 
They are almost not visible in the Case A simulation, with the 
wind of r] A in the — x direction looking continuous, but containing 
shells of enhanced density at periodic intervals. Fig. B5 also shows 
that the gradual fragmentation of the shell of dense r/ A material on 
the +x side of the system produces a series of small, cold, dense 
blobs within the larger, low-density wind cavity created by r) B . The 
higher the value of M Va , the larger and denser these blobs are, and 
the longer they last before mixing with the surrounding ^ B wind. 
Fig. B6 illustrates how in Case A, the winds on the — x side of the 
system in the xz plane remain cold, while the wind cavity on the +x 
side remains hot. In Cases B and C, the low-density cavities on both 
sides of the system remain hot, with only the intermediate regions 
of dense r/ A wind cold. 

In Fig. B7, we see that the extended narrow wind cavities on 
the +y side of the system in the yz plane are always smaller than 
the wind cavities carved on the — y side. These smaller cavities 
are always cold, while the larger cavities are always hot (Fig. B8). 
The difference in size between these cavities grows with decreasing 
M j]a . In Case A, the cavities are of comparable size, but in Case C, 
the cavities on the — y side are much larger. The temperature of the 
gas in the —y cavities is also nearly an order of magnitude hotter 
in Cases B and C than in Case A. Finally, the shells of dense r] A 
material on the —y side of the system are thicker and remain intact 
longer in the higher M, 1A simulations. 
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Figure Bl. Same as Fig. 1 . but for slices in the xz plane containing the binary semimajor (x) and orbital angular momentum (z) axes. 
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Figure B2. Same as Fig. Bl, but with colour showing log temperature. 
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Figure B3. Same as Fig. 1, but for slices in the yz plane containing the binary semiminor (y) and orbital angular momentum (z) axes. 
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Figure B4. Same as Fig. B3, but with colour showing log temperature. 
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Figure B5. Same as Fig. 8, but for slices in the xz plane containing the binary semimajor (x) and orbital angular momentum (z) axes. 
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Figure B6. Same as Fig. B5, but with colour showing log temperature. 
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Figure B7. Same as Fig. 8, but for slices in the yz plane containing the binary semiminor ( y ) and orbital angular momentum (z) axes. 
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Figure B8. Same as Fig. B7, but with colour showing log temperature. 
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